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Q^- ABSTRACT 

o : 

, The existence of a weak galactic magnetic field has been repeatedly confirmed by observa- 

^ ' tional data. The origin of this field has not as yet been explained in a fully satisfactory way and 

• • , represents one of the main challenges of the astrophysical dynamo theory. In both the galactic 

. ^ I dynamo theory and the primordial-origin theory, a major influence is exerted by the small-scale 

' magnetic fluctuations. This article is devoted to constructing a systematic second-order statis- 

tical theory of such small-scale flelds. The statistics of these fields are studied in the kinematic 
approximation and for the case of large Prandtl numbers, which is relevant for the galactic and 
protogalactic plasma. The advecting velocity field is assumed to be Gaussian and short-time cor- 
related. Theoretical understanding of this kinematic dynamo model is a necessary prerequisite 
for any prospective nonlinear dynamo theory. The theory is developed for an arbitrary degree 
of compressibility and formally in d dimensions, which generalizes the previously known results, 
elicits the structure of the solutions, and uncovers a number of new effects. The magnetic energy 
spectra are studied as they grow and spread over scales during the initial stage of the field ampli- 
fication. Exact Green's-function solutions are obtained. The spectral theory is supplemented by 
the study of magnetic-field correlation functions in the configuration space, where the dynamo 
problem can be mapped onto a particular one-dimensional quantum-mechanical problem. The 
latter approach is most suitable for the description of the kinematic dynamo in the long-time 
limit, i.e. when the magnetic excitation has spread over all scales present in the system. A simple 
way of calculating the growth rates of the magnetic fields in this long-time limit is proposed. 
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Subject headings: galaxies: magnetic fields — ISM: magnetic fields — magnetic fields — meth- 
ods: analytical — MHD — turbulence 

1. INTRODUCTION 

1.1. Astrophysical Motivation 

The question of the origin of the galactic magnetic field has long been a subject of much interest in 
plasma astrophysics. The existence of galactic magnetic fields was first inferred by Alfven (1937a, b) and 
Fermi (1949) from the properties of cosmic rays. This was later confirmed by observational data (Hiltner 
1949; Hall and Mikesell 1949). Modern observations indicate that our Galaxy possesses a magnetic field that 
has a large-scale component that is several /xG strong and is coherent on the scales of up to a kiloparsec 
(^ 10^^ cm). These scales are intermediate between the diameter of the galactic disk (~ 10 kpc) and its 
thickness (~ 100 pc). Observations indicate that magnetic fields of similar magnitudes (which correspond 
to magnetic-energy densities comparable to those of the fiuid motions of the interstellar matter) and spatial 
coherence are common in other galaxies as well (see reviews by Kronberg 1994; Beck et al. 1996; Zweibel 
and Heiles 1997). Most of the theories that have been advanced to explain the presence of these fields have 
in one way or another connected their origin with the dynamo action of the interstellar turbulence. 

The interstellar medium (ISM) consists of a partially ionized plasma that is regularly stirred on the 
scale of about 100 parsec by the shock waves generated by supernova explosions. The Reynolds number of 
the ISM is of the order of 10^, which allows a fully developed Kolmogorov-type turbulent cascade to be set 
up. The energy generated at the outer scale of 100 parsec is thus transfered approximately four decades 
down to the Kolmogorov inner scale, where it is dissipated by the molecular viscosity. If an initial seed 
magnetic field is introduced into such a medium, the turbulent velocity field should be expected to stretch 
the magnetic-field lines and thus amplify the field via the usual fast-dynamo mechanism (Sakharov 1982; 
Vainshtein and Zeldovich 1972). The turbulent nature of the physical processes involved clearly necessitates 
a statistical description. 

Traditionally, the hopes for a theoretical explanation of the galactic magnetic field have focused upon 
the mean-field dynamo theories (Parker 1955; Braginskii 1965b; Steenbeck, Krause, and Radler 1966; Moffatt 
1978; Parker 1979; Ruzmaikin, Shukurov, and SokoloflF 1988). The essential idea is to follow the evolution of 
the volume-averaged (i.e., effectively, large-scale) magnetic field subject to two main assumptions. First, it is 
assumed that the initial seed field is small and thus, during the initial stage of its evolution, the magnetic field 
is too weak to exert a significant amount of back reaction on the hydrodynamic motions. The Lorentz forces 
can therefore be neglected and magnetic field considered passive, which gives rise to the so-called kinematic 
approximation. Second, the large-scale mean field is assumed to be much stronger than the small-scale 
magnetic fluctuations.^ It turns out that, if these two assumptions hold, and if the interstellar turbulence 
lacks mirror invariance, the mean field will grow exponentially at a rate proportional to the amount of helicity 
possessed by the turbulent medium. Such mean-field amplification has come to be referred to as the a effect. 

It must, however, be appreciated that the mean-field theory faces serious challenges to its theoretical 
and physical foundations. It was pointed out already by Batchelor (1950) that the initial amplification of 



^This assumption can be circumvented for the Gaussian (5-correlated velocity field discussed in § 1.2 (Vainshtein 1970; Kulsrud 
and Anderson 1992; Boldyrev 2001). 
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the magnetic field in the kinematic regime would be accompanied by the transfer of the magnetic energy to 
small (nonhydrodynamic) scales. A wide range of such scales exists as a result of the huge disparity between 
the hydrodynamic and the magnetic Reynolds numbers of the ISM: while Re ~ 10^, its magnetic counterpart 
can be as high as Km ^ 10^^. The ratio of these two numbers, commonly referred to as the magnetic Prandtl 
number, Pr ^ R„i/Rc, determines the width of the scale interval between the Kolmogorov inner scale fc^, 
where the viscous dissipation cuts off the hydrodynamic turbulence spectrum, and the resistive scale krj, 
where magnetic energy is dissipated by the Ohmic resistive-diffusive mechanism. Since kj^/ki, ~ Pr^^^, 
the subviscous scales accessible to the small-scale fluctuating magnetic fields extend over as many as seven 
decades in the wave-number space. In the kinematic approximation, the rates of the exponential growth of 
the small-scale magnetic-fluctuation energy and of its transfer toward the small scales turn out to greatly 
(by the factor of 10^) exceed the growth rate of the mean field due to the a-effect, which operates on a 
time scale associated with the overall galactic rotation. Unlike the mean field, the magnetic energy grows 
regardless and independently of whether the turbulent velocity field has a helical component. Thus, the 
energy of the small-scale magnetic fluctuations must be expected to grow to equipartition with that of the 
smallest eddies of the turbulent velocity field much faster than the mean field can reach any appreciable 
values. The validity of the kinematic mean-field theory will therefore break down long before the growth 
mechanism predicted by it has time to manifest itself (Kulsrud and Anderson 1992). 

The mean-field approach is an attempt to explain the large-scale galactic magnetic field in terms of a 
coherent volume-averaged field being amplified by the hydrodynamic helicity effect. Besides the impossibility 
to neglect the small-scale fields, there exists another physical consideration that makes justifying this view 
rather problematic. Since no magnetic-energy dissipation mechanism is available at the large scales, the 
growth of the large-scale field must be consistent with the fiux-conservation constraint. The most popular 
theory has been that the magnetic-field lines are partly expelled from the galactic disk (see Parker 1979; 
Kulsrud 1999, 2000). However, the possibility of flux expulsion has been increasingly in doubt (Raflkov and 
Kulsrud 2000). 

In view of the fundamental difficulties associated with the mean-field approach, one is hard-pressed to 
seek alternative ways to construct an adequate galactic dynamo theory. A promising avenue of investigation 
naturally presents itself in the context of the small-scale-field amplification. It seems quite reasonably clear 
and has been an accepted point of view since the work of Batchelor (1950) that the growth of the fluctuating 
flelds at subviscous scales should culminate in the magnetic-fluctuation energy equalizing with that of the 
smallest turbulent eddies. What happens next is a subject of much interest and disagreement. The basic 
point of contention is whether the magnetic energy would saturate at viscous (Batchelor 1950) or resistive 
(Vainshtein and Cattaneo 1992; Gruzinov and Diamond 1994, 1995, 1996) scales, or rather proceed to reach 
full equipartition with the turbulence (Biermann and Schliiter 1951). 

An explanation of the existence of the large-scale galactic magnetic field could be within grasp if it were 
to be demonstrated that, once the magnetic-fluctuation energy had equalized with that of the smallest eddies, 
an inverse cascade would take place, forcing the magnetic excitation to keep growing and to gradually move 
towards ever-larger scales until equipartition were achieved at all scales up to the energy-containing ones. A 
statistical theory of this sort would predict the emergence of large-scale magnetic fluctuations with energy 
comparable to that of the turbulence (see Sec. 4 for further discussion of these matters). From the point 
of view of the mean-field approach, these fluctuations would still have a zero statistical average. However, 



^In this context, it should be noted that resolving such a broad range of scales in a numerical simulation is all but impossible, 
so theoretical understanding of the physics of the small-scale fields is indispensable. 
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observationally, they would locally appear as a field coherent at large scales (cf. Blackman 1998). 

This idea faces a serious setback when measured against the actual parameters of the ISM. Namely, the 
typical energy-containing scales of the galactic turbulence, to which one can expect the inverse cascade to 
bring the magnetic energy, are of the order of 100 parsec (the average distance between the supernovae, which 
drive the turbulence), or about 10 times smaller than the coherence scale of the observed large-scale galactic 
magnetic field. This problem can be resolved if one follows Kulsrud et al. (1997a) in their recent suggestion to 
consider the possibility of a primordial origin of the galactic magnetic field. The primordial-origin hypothesis 
regards the large-scale galactic field not as a product of galactic turbulent dynamo mechanisms, but as a 
residue of the analogous processes that take place in protogalaxies before they collapse into galaxies. 

Indeed, the conditions in the protogalactic plasma cloud seem to be more favorable than those in the 
galactic ISM for the operation of a bona fide turbulent dynamo. As does the ISM, the protogalactic plasma 
supports a Kolmogorov-type turbulence, this time driven by the shock waves originating from the instabilities 
associated with the forces of gravitational collapse. These instabilities occur on the scales comparable to the 
size of the protogalaxy (^ 100 kiloparsec) and thus, unlike galaxies, protogalaxies know no disparity between 
the system size and the energy-containing scale of the turbulence. The hydrodynamic Reynolds number of 
the protogalactic plasma is Re ~ 10"^, so an inertial range about three decades wide is available for the 
Kolmogorov turbulent energy cascade. Numerical simulations indicate that such a cascade is indeed set up 
(Kulsrud et al. 1997a). The protogalactic Prandtl number is even larger than the galactic one: Pr 10^^, 
which allows as many as eleven decades of subviscous scales accessible to the magnetic fields. 

The seed magnetic fields in the protogalaxy could be created from an initial state with no fields at all 
by the so-called Biermann-battery mechanism (Biermann 1950; Kulsrud et al. 1997a). This battery action 
is associated with the thermoelectric term in the plasma Ohm's law. This term is non-zero provided the 
motions are nonbarotropic (Vp x Vp ^ 0). Since the violent beginnings of the protogalax;y may be assumed 

to allow nonbarotropic large-scale pressure fluctuations, it can be estimated that the battery is capable of 
generating large-scale magnetic fields of about 10^^^ G before the dynamo term becomes dominant. 

The magnetic fields are then amplified and carried over to small scales by the same Batchelor (1950) 
mechanism that foiled the galactic mean-field theory. After the equipartition with the smallest turbulent 
eddies is reached, an inverse cascade may be envisioned that brings the magnetic energy back to the large 
scales. As the protogalaxy collapses into the galactic disk, its magnetic field is compressed and becomes 
the initial magnetic field of the newly formed galaxy. If the pre-collapse protogalactic magnetic field is 
amplified by the protogalactic turbulence to sufficiently large values, it may well be enough to give rise, 
upon compression, to a galactic field with the strength and scale of coherence reasonably close to those 
observed (Howard and Kulsrud 1997; Kulsrud et al. 1997b; Kulsrud 1999, 2000). [A midway approach is to 
consider the residual field resulting from the protogalactic dynamo as the seed field for further operation of 
the galactic dynamo mechanisms such as, e.g, the a-fl dynamo (Parker 1979; Kulsrud 1999, 2000) (see also 
Blackman 1998).] While this scenario is certainly largely hypothetical, it offers a physical possibility that 
can hardly be neglected. 

In all of the above discussion, the fluctuating magnetic fields in general and the small-scale fields in 
particular have played a prominent role. Under additional assumptions about the turbulent velocity field, 
the kinematic stage of their evolution can be treated analytically and exactly. Any theoretical understanding 
of the processes that occur once the nonlinear stage is initiated, must rely on a solid statistical theory of 
the small-scale kinematic dynamo. In this paper, we give a detailed and systematic exposition of the small- 
scale kinematic dynamo theory, following an often-neglected Cartesian principle: "... je me persuadai . . . 
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que, pour toutes les opinions que j'avais regues jusques alors en ma creance, je ne pouvais mieux faire que 
d'entreprendre une bonne foi de les en oter, afin d'y en remettre par apres, ou d'autres meilleures, ou bien 
les memes, lorsque je les aurais ajustees au niveau de la raison" (Descartes 1637).* 

1.2. The Kinematic Dynamo Model 

The full-scale nonlinear theory of the astrophysical dynamos is, of course, a formidable challenge and, in 
practical terms, a degree more complicated than even the problem of fluid turbulence. One must therefore 
try to find simplified models in order to be able to extract any theoretical understanding of the magnetic 
dynamo problem. The kinematic stage of the magnetic-field growth, when the field is so weak that it exerts 
no appreciable back force on the medium, has traditionally been the focus of theoreticians' attention. Indeed, 
if the velocity field is decoupled from the magnetic field and can be fed into the problem in some given form, 
the induction equation becomes formally linear (in the magnetic field) and the task of solving it starts to 
appear within theory's reach. Thus, we consider the magnetic field evolving according to the induction 
equation (formally in d dimensions): 

dtB' = -vJ'Bi^ + uij^B^ - u%B' + r]AB\ (1) 

where u*(t, x) is the externally specified advecting velocity field, it'^ = du^/dx^, B^j, = dB^/dx^, and the 
Einstein summation over repeated indices is used throughout. Two basic avenues of research have been 
pursued in the literature on the kinematic dynamo problem. One of them is to study the magnetic-field am- 
plification in some prescribed deterministic, though possibly chaotic (i.e. stationary and stochastic in space), 
fluid flows and/or in specific geometries (sec, e.g., Moffatt 1978; Ruzmaikin, Shukurov, and Sokoloff 1988; 
Childress and Gilbert 1995). Another approach is to consider random (in time) velocity fields with pre- 
scribed statistics and try to determine the resulting statistics of the magnetic field (Kazantsev 1968; Kulsrud 
and Anderson 1992). As the dynamo models of the latter kind appear to better refiect the turbulent nature 
of the fluid motions that occur in the astrophysical applications of interest to us (ISM and protogalaxies: 
see §1.1), we restrict ourselves to such statistical line of research only. Additional assumptions of spatial 
homogeneity and isotropy of the system and statistical stationarity of the velocity field are usually made. 
Furthermore, we assume mirror-invariance, thus excluding the helicity effects, which are not expected to be 
felt at the time scales relevant to the small-scale- field physics (see discussion in § 1.1 above and in Kulsrud 
and Anderson 1992). 

Unfortunately, the kinematic assumption alone is not sufficient to turn the turbulent dynamo problem 
into a tractable one. The realistic turbulent velocity fields possess nontrivial intermittent statistics, which 
are hard to handle. It is necessary to assume something about them, fully appreciating the fact that any 
such assumption is bound to be highly artificial. An obvious such assumption, to which many turbulence 
theories gravitate, is to consider a Gaussian velocity field. This brings about an enormous simplification due 
to the splitting property of the Gaussian averages. However, even the Gaussian assumption does not fully 
remove the usual turbulence moment closure problem. The remaining complication has to do with the fact 
that, in general, the advecting flow is finite-time correlated. 

When the advecting velocity field has a finite correlation time, the equations for statistical quantities 



I became convinced . . . that the best I could do with all the notions that I had thitherto taken into my credence was 
to have the good sense of ridding my mind thereof, in order to subsequently put back in either other, better ones, or even the 
same, once I had arranged them according to reason" (Descartes, Discourse on Method). 
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such as the moments of the magnetic field cannot be obtained in a closed form. Moreover, introducing a 
finite velocity correlation time into the problem leads to sensitive dependence of the resulting statistics on 
the specific structure of the velocity correlations in time and space. [For more detailed discussion of these 
issues and further results with regard to the finite-correlation-time effects, we address the reader to papers 
by van Kampcn (1976); Boldyrcv (2000); Schckochihin and Kulsrud (2001), and references therein.] In order 
to obtain a solvable kinematic dynamo model, we must therefore take another drastic simplifying step. This 
consists in assuming that the velocity field, besides being Gaussian, also possesses the white-noise property, 
i.e., it is a random process 6-correlated in time: 

(«,'(*. x)w^'(t',x')) =(5(t-i')K*^'(x-x'). (2) 

Such a synthetic velocity field is sometimes called the Kraichnan ensemble, in deference to the author who 
first proposed it as a model for studying the passive-scalar problems such as the advection of temperature, or 
of concentration of an admixture (Kraichnan 1968, 1994). In application to the magnetic fields, this model 
was proposed independently by Kazantscv (1968), whose work, in fact, predates Kraichnan's. In physical 
terms, the Kazantsev-Kraichnan assumption amounts to the short-correlation-time approximation and is 
valid if the advecting velocity field is correlated at times much smaller than the characteristic time of the 
dynamo action. The latter is of the order of the inverse velocity gradient (the eddy-turnover time of the 
velocity field). 

When the white-noise assumption is put in force, the closure problem resolves itself. Thus, one of the 
main reasons for the Kazantsev-Kraichnan passive-advection paradigm being so attractive is that it is one 
of the few available "Ersatz-models" of turbulence, which hold promise of exact solvability. However, the 
existing theory is not complete, and the statistics of passive advection continue to generate considerable 
interest, for the problem of passive advection can serve as a vehicle for developing a set of analytical tools 
appropriate for attacking more realistic turbulence models. 

The objective of this paper will be to develop a theory that would allow us to calculate, under the 
assumptions explained above, the two-point second-order correlation function of the advected field S*: 

{B\t,x)B^ {t,x')) = {t,x - x'). (3) 

Here and everywhere in this work, the angle brackets denote ensemble averaging with respect to the random 
advecting velocity field u% as well as with respect to the initial distribution of B\ The latter, however, 
should not greatly affect matters after a finite transient time has elapsed (physically, about one eddy- 
turnover time). The two-point, second-order correlation functions have a clear physical interpretation: their 
Fourier transforms are the energy spectra. 

1.3. General Scenario of the Small-Scale-Field Dynamo 

The study of the kinematic dynamo problem in the above formulation was pioneered by Kazantsev 
(1968). Kazantsev reduced the problem of finding the two-point correlation functions of the fluctuating 
magnetic fields to solving a certain quantum mechanical problem that described a particle with variable 
(position-dependent) mass moving in a one-dimensional potential well. The eigenfunctions and the energy 
levels of the particle's Hamiltonian corresponded to the magnetic-field correlation functions with stationary 
spatial profiles and their overall growth rates, respectively. The ground state of the particle gave rise to the 
largest growth rate and thus determined the long-time evolution of the magnetic fluctuations. Importantly, 
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the exact shape of the potential depended on the particular form of the velocity correlation functions.^ 

In order to extract the general properties of the problem, Kazantsev concentrated his attention on the 
fluctuating magnetic fields at scales very small compared to the correlation length of the velocity field. In 
this limit, the two-point velocity correlation functions, which depend on the distance between the points, 
can be expanded around the origin. Two principal subclasses of the two-point kinematic-dynamo problem 
can be identified in this context, depending on the physical setting one is interested in. If a problem with 
a large magnetic Prandtl number is considered, i.e., if the resistive-diffusive scale k^^ is much smaller than 
the characteristic scale of the advecting flow, viz., we find ourselves dealing with 

what is customarily referred to in the literature as the Batchelor regime [after the original work of Batchelor 
(1959) where this view of passive advection was first elaborated]. In this case, the Taylor expansion of the 
velocity correlator K*-'(y) can be used, whose lowest-order y-dependent term is quadratic (see Appendix A). 
Since in the astrophysical applications that are of most interest to us (the interstellar medium and the 
protogalactic plasma), the Prandtl number is very large, it is the Batchelor regime that we shall concentrate 
on in this work. The models with small Prandtl numbers, where the scale ratio between the flow and the 
magnetic field is reversed, kn/ki, ~ Pr^^^ < 1, lead to nonanalytic y dependence for the velocity correlator 
that magnetic fluctuations "feel" at small scales. Slightly different quantum mechanics result (Kazantsev 
1968; Vainshtein 1982a,b; Kichatinov 1985; Vainshtein and Kichatinov 1986; Vergassola 1996; Rogachevskii 
and Kleeorin 1997). These matters fall beyond the scope of this paper and will be reported elsewhere. 

Let us now outline the general scenario that emerges from the investigation of the statistical evolution of 
small-scale fluctuating magnetic fields in the Batchelor regime (Kulsrud and Anderson 1992). Let us assume 
that, initially, the magnetic fields are excited at scales comparable to the characteristic scale of the advecting 
flow k~^. (In the astrophysical context discussed above, this assumption is consistent with the current 
undertanding of the physics of seed fields — see § 1.1 and references therein.) Due to the convective and 
line-stretching ( "dynamo" ) action of the velocity field, two processes are initiated and proceed exponentially 
fast in time, at the rates comparable to the eddy-turnover rate of the velocity field. First, the magnetic 
fluctuations spread over a widening range of scales, their bulk shifting towards ever-smaller onc;s. A power- 
like spectrum forms behind the advancing peak. Second, the level of excitation (amplitude) of each spatial 
Fourier mode either grows (in 3D) or decreases (in 2D). The combination of these two processes results 
in the exponential growth of the total magnetic energy. Fig. 1 and Fig. 2 illustrate the evolution of the 
magnetic-fluctuation spectra in two and three dimensions respectively. Note that in 2D, the overall growth 
of the energy is due entirely to the spreading of the excitation over an increasing number of modes, while 
in 3D, each mode also grows by itself. We call this stage of the kinematic dynamo the diffusion-free regime. 

Eventually, the magnetic fluctuations reach the diffusive (resistive) scales ~ k~^. The effect of the 
Ohmic dissipation is to check the decrease of the characteristic scale of the magnetic field. The evolution of 
the magnetic energy is now determined exclusively by the exponential growth (in 3D) or decay (in 2D) of 
the level of excitation in each particular mode. 

We will see in our analysis that, while the evolution of the magnetic fluctuations starts at scales compa- 
rable to those of the advecting velocity (henceforth referred to as large, or integral, scales), practically all of 
their energy will shift to the small scales, k^ k„, after a few eddy- turnover times. Thereupon, the spectral 



^Recently it was demonstrated that statistics of general passive tensor fields, such as gradients of temperature or concen- 
tration, tensor products of vectors, etc., can also be reduced to a quantum mechanical problem. Remarkably, this problem is 
exactly solvable and describes d particles on the line which interact pairwisc with the potential ~ — xf)^, the so called 

Calogero-Sutherland potential (Bernard, Gaw§dzki, and Kupiainen 1998; Boldyrev and Schekochihin 2000). 
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properties of the small-scale fields can for the most part be understood in disregard of the small fraction of 
the magnetic energy that is situated at the large scales. The auxiliary model where the effect of the large 
scales is neglected and all the magnetic energy is considered to be concentrated at small scales, will also be 
called the diffusion-free regime before the resistive scales are reached, and the resistive regime afterwards. 
For both cases, exact Green's-function solutions can be obtained (see §2.3 and §2.4). 

The eventual spreading of the excitation over all scales available in the system, including the resistive 
and the large, ushers the true long-time asymptotic regime (or the long-time limit), which was the principal 
subject of Kazantsev's investigation. In three dimensions, a stationary correlation profile is expected to form 
which, if the Prandtl number exceeds a certain critical value, continues to grow exponentially. This is the 
cigcnstatc arising from Kazantsev's quantum mechanics, and Kazantsev's objective was to determine the 
growth rate and the scaling of the correlation function. 

While the analysis of the small-scale asymptotics allowed him to constrain this growth rate between 
certain, relatively narrow, bounds, Kazantsev identified the main difficulty of the theory: effectively, the 
potential that appeared in his quantum mechanics was of the inverse-square kind (^ —l/x^) and therefore 
had to be regularized at both small (resistive) and large (velocity-field) scales, in order for the cigcnstatcs 
of the problem to be determined (see, e.g., Morse and Feshbach 1953). While on the small-scale side the 
resistive regularization could be included exactly, the boundary conditions at large scales were unknown, 
so there was not enough information to fix the dynamo growth rate. A number of authors took up the 
study of Kazantsev's quantum-mechanical model in the 1980's. Much mathematical insight was gained into 
the properties of Kazantsev's quantum mechanics, but the integral-scales problem was mostly circumvented 
by choosing particular reasonable forms of the velocity correlation function (Ruzmaikin and Sokoloff 1981; 
Novikov, Ruzmaikin, and Sokoloff 1983; Artamonova and Sokoloff 1986; Kleeorin, Ruzmaikin, and Sokoloff 
1986; Maslova and Ruzmaikin 1987). 

The particular cases analyzed in the literature indicate that the second-order statistics (energy spectra 
and growth rates) are insensitive to the exact form of the large-scale regularization (see Novikov, Ruzmaikin, 
and Sokoloff 1983; Artamonova and Sokoloff 1986; Kleeorin, Ruzmaikin, and Sokoloff 1986; Maslova and 
Ruzmaikin 1987; Kulsrud and Anderson 1992; Subramanian 1997, 1998, 1999, §3.3, and Appendix C 
below). The growth rate in three dimensions can be determined exactly and turns out to be slightly smaller 
[by ~ 1/ (InPr)^] than the growth rate of each individual mode that operated in the diffusion- free regime. 
In two dimensions, the steadily growing bound eigenstates do not exist, and a continuous spectrum of 
eigenvalues is formed that is bounded only by zero and therefore indicates that the exponential decay of 
the magnetic fluctuations in the resistive regime is replaced by a power-like decay in the long-time limit 
(Novikov, Ruzmaikin, and Sokoloff 1983). In this context, we recall that, as was predicted most generally 
by Zeldovich's antidynamo theorem (Zeldovich 1957), only transitory growth is possible in two dimensions, 
and the magnetic field must eventually be dissipated by the resistivity. 

The goal of this paper is to give a coherent general exposition of the two-point second-order statistical 
theory of the kinematic dynamo. The emphasis is naturally placed upon the universal aspects of the problem, 
which emerge when the statistics of small-scale fields are considered. As long as specific functional features 
of the velocity correlation function (2) remain inessential, only its small-scale expansion affects the statistics 
of the passive magnetic fields: 

1^'' (y) = Ko5'^ -\k2 {y^5'i + 2ayY) +■■■ ■ (4) 

In this case, the kinematic dynamo problem contains only two essential dimensionless parameters: the di- 
mension of space d and the degree of compressibility a of the velocity field. Rather than considering various 
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particular cases, such as incompressible, a = — 1), irrotational, a = 1, three-, or two-dimensional 

dynamo, we work in general, explicitly keeping all dependences on these parameters. This enables us to 
construct a solid theoretical framework which incorporates all previously available results, elicits their un- 
derlying structure, and allows us to establish several previously unknown features and facets of the kinematic 
dynamo problem. Besides the general evolution equations for the magnetic-field spectra and configuration- 
space correlation functions that arc derived for arbitrary compressibility and spatial dimension [such as 
equations (11), (21), (56), and (61)], the new results include, for example, the exact Green's-function solu- 
tion for the small-scale-field spectra in the resistive regime [Eq. (35)] and a simple regularization procedure 
that allows to calculate the dynamo growth rates in the long-time limit (§ 3.3). 

This paper contains three more Sections and three Appendices. In Sec. 2, the kinematic dynamo problem 
is studied in its spectral form. Time evolution of the magnetic energy and its propagation over scales are 
described up to the point when all available scales are excited: from those where the velocity field operates 
to those at which the magnetic energy is diffused by the Ohmic resistive damping. In Sec. 3, we study 
the two-point correlation functions of the magnetic fluctuations in the spirit of Kazantscv's configuration- 
space quantum-mechanical approach, which is convenient for handling the kinematic dynamo problem at all 
scales in the long-time limit. Each of these sections starts with a brief summary of its internal organization 
and of the results contained therein. Sec. 4 contains concluding remarks of physical nature. In particular, 
we discuss the possible directions of future research and some strategies whereby the yet-to-bc-devclopcd 
nonlinear small-scale-dynamo theory could be built upon the foundation laid by the kinematic theory worked 
out in this paper. Under Appendix A, we include a number of technical appendices dealing with the basic 
properties of the correlation functions of isotropic fields. In Appendix B, we state the fundamental Gaussian 
averaging theorem that lies at the core of the derivations of closed equations for the correlation functions 
in the Kazantsev-Kraichnan model: the Furutsu-Novikov formula (Furutsu 1963; Novikov 1965). Finally, 
Appendix C contains some overflow asymptotic results pertaining to the material of Sec. 3. 



2. SPECTRA OF MAGNETIC FLUCTUATIONS 

We will start our investigation of the kinematic dynamo by looking at the spectra of magnetic fluctua- 
tions. While in most of the works on the subject (Kazantsev 1968; Ruzmaikin and Sokoloff 1981; Molchanov, 
Ruzmaikin, and Sokoloff 1983; Novikov, Ruzmaikin, and Sokoloff 1983; Kichatinov 1985; Vainshtein and 
Kichatinov 1986; Artamonova and Sokoloff 1986; Klccorin, Ruzmaikin, and Sokoloff 1986; Maslova and 
Ruzmaikin 1987; Vergassola 1996; Rogachevskii and Kleeorin 1997; Subramanian 1997, 1998, 1999), the 
configuration-space view was chosen because of a transparent quantum-mechanical form the problem took, 
it is easier to gain a general physical understanding of the processes under way if one first gains a firm 
foothold in the wave-number space. The statistical-evolution scenario outlined in the Introduction, viz., the 
growth and spreading of the magnetic excitation over scales, will naturally emerge in the k-space description. 
Another important consideration is that the correlation functions are better recognized as such in the Fourier 
space, while in the x space, only the functions that are Fourier transforms of the proper k-space correlation 
functions constitute the class of allowed solutions (see Appendix A.l). 

The evolution equation for the magnetic-field spectrum was first derived by Kazantsev (1968) and, in 
a somewhat different, but equivalent, form, by Kraichnan and Nagarajan (1967). Kulsrud and Anderson 
(1992) later developed a rather detailed theory of magnetic fluctuation spectra. Here we derive some of the 
existing results in a more general form (in d dimensions and for arbitrary degree of compressibility of the 
advecting flow) and then proceed to present several new results on the spectral theory of small-scale fields. 
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in § 2.1, we derive the general integral equation that governs the evolution of the magnetic spectrum. From 
it, we infer the exponential growth of magnetic energy and the exponential decrease of the characteristic 
scale of magnetic fluctuations. In § 2.2, we obtain a partial differential equation that describes the behavior 
of the magnetic spectra at scales much smaller than those of the advecting flow. In § 2.3, we study the 
exact Green's-function solution for the spectrum of the small-scale magnetic fluctuations in the diffusion-free 
regime. In § 2.4, we proceed to the resistive regime. The previously unknown exact Green's-function solution 
is obtained in quadratures that is valid simultaneously in both inertial and resistive scale ranges. In § 2.5, we 
take up the question of the long-time asymptotics of the magnetic spectra that emerge after both resistive 
and integral scales have been excited. The latter problem is better suited to be treated in the configuration 
space. This brings us to the material presented in Sec. 3. 



2.1. The General Equation for the Magnetic-Energy Spectrum 



We shall consider the kinematic dynamo problem as a problem of Gaussian passive advection defined 
by equations (1) and (2). The velocity correlation function in the wave-number space is, by definition. 



(5) 



(in the case of incompressible flow, k = —k, while for the purely irrotational flow, k = 0). We will also 
sometimes use the spectral energy function I{k) = K"{k) = dK{k) + k{k). The magnetic correlation function 
is deflned as follows: 



{B\t,k)B^ {t,k')) = {2TrfS{k + k')W^{t,k), 



H'^{t,k) 



(6) 



The induction equation (1) in the Fourier space can be written as follows (from here on, we suppress the 
time dependence in the arguments unless it is essential for the clarity of exposition): 



dtB\k)+f^k^B\k)^tki J ^ [u'{k')BHk-k') - u\k')B\k-k')]. 

We can now write an equation for the tensor B'^(t,k)B^t,k') and average it: 

dt {B'{k)Bi{k')) + 7? (fc2 + k'^) {B\k)B^{k')) = 

iki 



(7) 



f d'^k" \ 
J '{2nY\ 



ik', 



{u\k")B^{k - k")B^{k')) - {u^{k")B\k - k")B^{k')) 

{u\k")B'(k' - k")B'{k)) - {u^{k")B\k' - k")S'(k)) 



(8) 



The mixed triple averages that have arisen in the right-hand side are split by the Furutsu-Novikov (or 
"Gaussian-integration" ) formula (see Appendix B) into products of velocity correlation functions defined in 
the formula (5) and averaged response functions: 

{B^{tM)B^{tMW{t,ks)) = f'it'J<i'k'{u\tM)u-{t\k'))(^^l^^^^^ 

' 5[Bi{t,ki)B^{tM)]\ 



- (27r)''«'"(k3) 



5u^{t, -ka) 



(9) 
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The averaged same-time response function that appears above is 

/5[B'(ki)BJ(k2)]' 



(5u"(-k3) 



6'"H^^ (ki + ka) - ^""if '^ (ki + ka) (5(ki + k2 + kg) 



2m 



^jn^mi(k2 + ka) - (5""if'^(k2 + ka) (5(ki + k2 + ka), 



(10) 



where the tensor (k) is defined by the formula (6). The expression (10) is obtained by formally integrating 
the (unaveraged) equation for B^{t,ki)B^ {t,k2) up to time t, taking the functional derivative 6/Su"'{t',]s.s), 
averaging, setting t' = t, and using the causal property of the response functionals. 

Using the definitions (5) and (6) and after some straightforward algebra, we arrive at the so called 
mode- coupling equation: 



d-1 



dui 



(27r)" 



dtH + 2{r]T + rj)k''H = 

'(c^-i)feVn+^('-'-')(^ '-").(n 



+ 



■,2 (k-k')^^ 



id -3) 



(k ■ 



K{k") + (d - 3) 



[K{k") + k{k")\ 
(k-k")2 



k'" 



k{k") \H{k'), 



(11) 



where k' = k — k", and we have introduced an auxiliary quantity that one may term the turbulent diffusivity: 

VT = -J-^,I{k) (12) 



[recall that I{k) = dK{k) + K{k)]. Eq. (11) is the generalization of the mode-coupling equations of Kazantsev 

(1968), Vainshtcin (1982a), and of Kulsrud and Anderson (1992) to the case of d dimensions and arbitrary 
degree of compressibility. It is now straightforward to verify that Eq. (11) integrated over all wave numbers 
gives the following evolution law of the total magnetic energy W = {2tt)~'^ J d'^kH{k): 



dtW = 2jW - 2r] k^ W, 



where the growth rate is 



27: 



d-l 
d 



d'^k 
{2nY 



k'^[2K{k) + k{k)] 



and the mean square wave number of the magnetic spectrum is 

d-^fe 



(13) 



(14) 



(15) 



The non-conservation of the magnetic energy reflected by Eq. (13) is understandable, because the system is 
not closed and the passive magnetic field constantly receives energy from the velocity field. Note that the 
white-in-time velocity field formally constitutes an infinite reservoir of energy. 

It is instructive to have also an evolution equation for fc^. This is derived in quite the same fashion as 
Eq. (13) and reads 



dtk^ = D + 272^2 - 2?7 k^- {k'^f 



(16) 
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where the growth rate is 



1 f d'^k 

272 = 



I{k)+4{K{k) + p-lk{k) 



(17) 



d J {2TTy 

the source term, which corresponds to a slow diffusion of the magnetic spectrum in k space, is 

^= V/(£^^'[''^(^^+^^^)]' ('^^ 

and the "fourth moment" of the magnetic spectrum is defined analogously to the second (15): 

All of the velocity-spectrum moments introduced above [formulas (14). (17), and (18)] can be shown to 
be positive based on the realizability properties of the isotropic velocity-field correlation functions (see 
Appendix A. 3). Note that, since k^ — (fc^)^ = {k^ — fc^)^ > 0, the resistive term still plays the ultimately 
balancing role in Eq. (16). 

Eq. (13) and Eq. (16) provide us with a basic understanding of the time-evolution properties of the 

magnetic energy spectrum. Let us assiimc that initially the magnetic fluctuations are concentrated at wave 
numbers comparable with or greater than the characteristic wave number ki, of the velocity field, and that 
the Prandtl number is large, i.e., the characteristic wave number associated with Ohmic resistive diffusion, 
kjf ^ Vv^^^kv, greatly exceeds k^. We see that the mean square wave number of the magnetic fiuctuations 
increases exponentially in time with the growth rate 72, which is of the same order as that of the magnetic 
energy, until this increase is checked by the resistive dissipation [Eq. (16)]. Thus, the bulk of the magnetic 
spectrum will shift towards the resistive scale range exponentially fast with characteristic rate 72. The 
fluctuations themselves will also increase in strength, exponentially fast as well, with the rate 7 [Eq. (13)]. 

Naturally, the question arises, what is going to happen when the resistive scales are finally reached. That 
is where the principal technical difficulties of this theory lie. Before we attack this problem, let us probe 
a little further into the additional information Eq. (11) can yield regarding the properties of the spectra of 
small-scale fields. 



2.2. The Small-Scale-Field Equation 



We saw in the previous section that the spectrum of magnetic fiuctuations had a tendency to shift to 
scales much smaller than those of the velocity field. This allows us to study the magnetic energy spectra 
in what we shall call the small-scale approximation. Namely, we may notice that the kernel of the integral 
operator in the right-hand side of the mode-coupling equation (11) essentially cuts off the range of values 
that k" may take at the scales k" ~ k^, and hence, for k^ k^.we may expand 



H{k') = H{k) + 



1 fk"y .2^ k" 

- — sm C — cos( 

2 \ k k 



kH'{k) 



+ 



cos" eeH"{k) + 0[{k"/kf 



(20) 



1 

2\k 

where 6 is the angle between k" and k. Expanding similarly the kernel and performing the necessary angle 
integrations (with the aid of Appendix A. 5), we arrive at the following small-scale- field (SSF) equation: 

dtH + 2'nk'^H = 7 {Ak'^H" + BkH' + CH) , (21) 
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where 7 is defined by the formula (14), while the dimensionless coefficients A, B, and C depend only on the 
dimension of space and the degree of compressibility of the velocity field: 

l + 2a ^^^^ 



(d- l)(da + 2)' 



_ d-l + Aa 

~ (d-l)(da + 2)' ^^"^^ 
^ 2(rf-2) 

(d-l)(da + 2)' ^ ' 

Here a is the so-called compressibility parameter of the velocity field which arises naturally the small- 
scale expansion of the velocity correlation function and is ubiquitous in one-point statistical theories (see 
Appendix A. 4 and Boldyrev and Schekochihin 2000; Schekochihin and Kulsrud 2001): 



1 + 



(25) 



By using the realizability properties of the velocity correlation functions (see Appendix A. 3), we find that 
the values of a are sandwiched between —l/{d+l) (incompressible flow) and unity (irrotational flow), A > 0, 
the sign of B is undetermined, and C > except in 2D, where it vanishes. The values of A, B, and C in the 
four special cases that are usually studied in the literature are given in Table 2. In the case of incompressible 
flow, Eq. (21) was earlier obtained in 3D by Kazantsev (1968) (see also Vainshtein 1980, 1982b; Kulsrud and 
Anderson 1992) and in d dimensions by Gruzinov, Cowley, and Sudan (1996). 

It is worthwhile noticing that A{d+l)d + Bd + C = 2, and therefore the SSF equation (21) satisfles the 
exact energy evolution law (13). Thus, the small-scale approximation is conservative. This is an essential 
property because it means that all of the magnetic energy is concentrated at the small scales. It is therefore 
reasonable to expect that the spectral properties of the small-scale fields would mostly be captured by the 
SSF equation (21). This consideration encourages us to embark on a detailed investigation of the solutions 
of this equation. 



2.3. Spreading of Magnetic Fluctuations in the Wave-Number Space: Diffusion-Free Regime 

First, let us pick a wave number fco such that <^ ko <^ kr,, and inquire how an initial spectrum 
concentrated at that wave number will evolve over time. This is easily determined, because neglecting the 

diffusive term in Eq. (21) and changing to the logarithmic variable z = Ink, renders this equation a one- 
dimensional heat equation in a moving frame. The Green's function for such an equation is well known, so 
we find that an initial spectrum such that Ho{k) oc 6{k — ko) spreads out over time into a widening lognormal 
profile: 

M{t,k) = -^k''-'H{t,k) 

Wo^XojtfkY" 1 [Hk/ko)r 
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where Sd = 2tt'^/'^ /r{d/2) is the area of the unit sphere in d dimensions, Wq is the initial magnetic energy, 
and we have introduced two new parameters: 

, A-B2-d~2a 
^0 = d-l + so, where = ^ = ^(i + 2a) ' ^'^^ 

2 (3 - d){d - 1 + 4a) + 4(rf - 2)(1 + 2a) - (1 + 2a)^ 
Ao = C-^So= 4(rf-l)(l + 2a)(rfa + 2) ^^^^ 

= 2-^(^0 + 1)'. 
The values of sq, ^o, and Aq in the standard cases are given in Table 2. 

Examining the solution (26), we immediately conclude that the following basic processes are initiated: 

(i) the fluctuation strength in each logarithmic wave-number band grows or decays exponentially at the 
rate A07, (ii) the number of such bands that are excited [i.e., the width of the lognormal envelope in the 
solution (26)] grows exponentially fast at the rate 4A7, (iii) a magnetic spectrum with exponent ^0 is formed 
inside the lognormal envelope, with the peak of the spectrum moving rightwards (to smaller scales/larger k): 

kpcak = ko exp {2£^o A-ft). It is not hard to verify that the effective sum of these processes produces the total 
magnetic-energy growth rate we have previously obtained [Eq. (13)]: 27 — A07 + A{^q + 1)^7. 

Let us make several observations regarding the behavior of the parameters of the solution (26) defined 
by formulas (27) and (28): the spectral exponent sq (or ^0), and the fraction Aq of the total growth rate 
that corresponds to the growth or decay of each individual Fourier mode. (The latter is plotted in Fig. 3 as 
a function of the compressibility parameter a in two and three dimensions.) 

(i) In 3D, the spectral slope that forms behind the advancing peak of the spectrum does not depend on the 
degree of compressibility: ^0 = 3/2 (sq = —1/2) regardless of the value of a. 

(ii) While the number of excited modes increases exponentially fast regardless of the particular parameter 
values, individual modes grow in 3D (Ao > for all a) and decay in 2D (Aq < for all a 7^ 0). 

(iii) In 2D and for a = 0, Aq = and sq = 0. Curiously, this suggests the forming of an equilibrium-like, 
equipartion spectrum {^q = d — 1). 



2.4. Spreading of Magnetic Fluctuations in the Wave-Number Space: Resistive Regime 

The solution (26) will cease to be valid when the magnetic excitation that has resulted from it reaches 
the integral scales ~ k„ and/or the resistive scales ~ fc^. Let us first consider the simpler of these two aspects 
of the problem, namely that where the resistive scales are reached before the integral ones. In this case, the 
further behavior of the magnetic fluctuations will be determined by the SSF equation (21) with the resistive 
term retained. 

When the diffusive term is allowed to reappear in the SSF equation (21), a Green's fimction solution 
can still be obtained in quadratures. In order to accomplish that, we seek the solution in the form 

H{t,k) = e^°'^'k'°f{t,k/kr,), (29) 

where Aq and sq are deflned by the formulas (28) and (27), and fc^ = (7^/277)^/^ is the characteristic wave 
number of the resistive dissipation. It is then elementary to see that f{t,x) satisfies 



dtf = A^{x^f" + xf'-x^f) 



(30) 
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The Green's function solution of this problem can be found by applying to Eq. (30) the Kontorovich-Lebedev 
transform (see, e.g., Erdelyi 1954), which, in the form adapted to our problem, can be defined as follows: 

/■°° da; 

F{t) = / ^KiAx)f{x) (direct), (31) 
Jo X 

f[x) = -o / dTTSvah.{'KT)Kir{x)F(T) (inverse), (32) 



where Kir{x) is the modified Bessel function of the second kind (also known as the Macdonald, or Basset, 
function). Kontorovich-Lebedev transforming Eq. (30), we get 

dtF{t,T) = -Ajr^F{t,T). (33) 

The solution of Eq. (30) is therefore 

2 f°° dr' f°° 2 

f{t,x) = — —hW) (iTTSmh{-KT)Kir{x)Kir{x')e-^'<^ \ (34) 
Jo X' Jo 

where fo{x) is the initial profile. To compare this result with the inertial-range solution (26), let us recast it 
in the form of an evolution law for a spectrum that is initially a 5-like spike, Mo(f, k) oc S{k — ko). We have 

^0 fk_Y" A r ........... ^. K fh\ .-A.rH 



It is not hard to verify that, for fco <C fc^ and fc <C fc^, the solution (26) is recovered as an asymptotic of the 
more general solution (35). On the opposite end of the spectrum. A: :» fc^, one finds the expected exponential 
cutoff due to resistive dissipation. Fig. 4 illustrates the evolution of the spectrum (35) with time. 

Integrating the spectrum (35), wc find that the total magnetic energy evolves according to 



W{t) = Wo e^"^' ( 



-(«o+l) 



2^0 r 
^Jo 



X / dTrsinh(7rT) 



r 



Co + 1 + 



K,r ('^') e-^-^^*. (36) 



This formula reflects the gradual slide of the "effective growth rate" of the magnetic energy, 7eff(i) = 
ln[VF(t)/VFo] , from the diffusion-free Hmit 7eff(i — > 0) = 27 to the diffusion-reduced value, 7eff(t — » 
00) = A07. This is demonstrated in Fig. 5. 

Direct integration of the solution (35) also produces an exact formula for k'^{t). Its asymptotic behavior 
can, however, be inferred already from that of the energy: 

F(t^oo) = ^^fc2^(Co + l)X (37) 

[see formula (28)]. This is a good quantitative measure of the characteristic wave number in the resistive 
range around which the bulk of the magnetic fluctuation energy will stabilize. The evolution of [fc^(f)] is 
plotted in Fig. 6. 

Thus, the solution (35) describes how the spreading lognormal spectrum hits the resistive scale range 
and and how any further refinement of the magnetic-fiuctuation scales is suppressed by resistive dissipation. 
The total magnetic energy now concentrates in the resistive scale range and grows (in 3D) or decays (in 2D) 
at the rate A07. This regime persists until the excitation finally reaches the integral scales. 
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2.5. Spectra and Growth Rates of Magnetic Fluctuations in the Long-Time Limit 

In order to gain some idea of the long-time behavior of the magnetic spectra, which sets in when all 
scales available in the system, including the integral ones, are excited, let us seek the solution of Eq. (21) in 
the form 

H{t,k) = e^'^*h{k). (38) 

The constant spectral profile h{k) must then satisfy the following linear ordinary differential equation of the 
Bessel kind: 

Ak^h" + Bkh'+(c-X-^k^^h = 0. (39) 



Before considering the explicit solutions of this equation, it is instructive to recast it in a form mathematically 



equivalent to that of a Schrodinger equation for a particle with energy equal to —1/k^ moving in a one- 



dimensional inverse-square potential, viz., 
where the "wave function" is 



^(fc) = k-'°+^/'^h{k), (41) 

So and Aq are defined by the formulas (27) and (28), and kr, = (7^/2//)^/^ is the inverse resistive scale. 
The potential in Eq. (40) is repelling and therefore cannot have a negative energy level for A > Aq + A/ 4. 
Solutions of the form (38) do not exist in this case. For A < Aq, the so-called "fall on center" occurs (see, 
e.g.. Landau and Lifschitz 1977), whereby the solution becomes oscillatory with infinitely many nodes as 
k approaches zero, so it cannot be a spectral function: for k <^ kn, 




^ In A + const). (42) 
kn J 

1/2 

Of course, the small-scale approximation is only valid for k ^ ki, A;,,Pr ' , so for A lying below Aq, but 
within some 5\ 0(1/ In^ Pr) of it, the solution will not yet have any nodes in the range of wave numbers 
for which Eq. (40) is valid. Thus, the growth rate of the stationary-profile solution (38) is contained in the 
interval 

Ao - ^A < A < Ao + A/4 = Amax, where ^ A ~ [ ■ (43) 

V In Pr / 



For the 3D incompressible case, the bounds (43) on the growth rate of a stationary-profile spectrum, 

as well as the solution (44), were first established by Kazantsev (1968). He also pointed to the fundamental 
problem that, since the boundary conditions at small k (integral scales) are unknown, there is no rigorous 
way to fix the value of A within the framework of the small-scale approximation. 

The two fundamental solutions of Eq. (39) are the modified Bessel functions Iv{x) (k/kr,) and K^(^x) {k/k^) 
times k^° , where v{X) = -^/(A — Xoj/A. The obvious requirement that the physically acceptable solution must 
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decay at large k, leaves us with 



M{t,k) 



Sd 



(2^) 



-^k^-^H{t,k)=C{X)e^'^^ 



(44) 



2«o-ir ( Co + l + KA) \ p + 1 - 



n -1 



where Ki, is the Macdonald function. The "subrcsistive" asymptotic of this sohition does not depend on A 
in any essential way and simply reflects the expected exponential cutoff of the magnetic spectrum; 



M(t,fc^oo)^y|c(A)e^^* (^^ 



?o-l/2 



k 

exp 1 - — 



The asymptotics valid for k <^ kr, are 
M{t, k) 



\kjj 



?o-i^(A) 



, A > Ao, 



M{t,k) ^ c{x)e^'',t(^k_y^^i\y ^^^^^ 



and, for A = Ao — 6X, where < ^A <C 1, we have 



M(..).;f.(A)e-(A)%J;f._^^ 



(45) 

(46) 
(47) 

(48) 



which is, of course, consistent with (42). 

Clearly, if the initial distribution of the magnetic fluctuations is chosen in such a way that it spreads 
over to resistive scales before it does to the integral ones (i.e., if the initial characteristic wave number ko in 
the solution (26) of §2.3 lies "closer" to fc^ than to fc^), there will be an intermediate period between the 
time the resistive scales are reached and the time when the excitation finally arrives at fc ~ fcj/. During this 
period, the boundary condition at the integral scales cannot, of course, affect the spectrum. The solution of 
the form (38) consistent with the spectra obtained in § 2.3 and § 2.4 [solutions (26) and (35)] wiU clearly be 
the one with A = Ao: indeed, it can be checked that the solution (44) with A = Aq is the long-time asymptotic 
of the solution (35). This is essentially what was demonstrated numerically by Kulsrud and Anderson (1992) 
(see also the recent work of Chertkov et al. 1999). Our Green's- function solution (35) is the mathematical 
expression of this fact. 

We may further argue that, if we were to mandate that the solution (44) vanish at some ko such that 
ki, <^ ko <^ kr,, then, by virtue of the asymptotic expression (48), we would be required to set 



A = Ao- 



[ln(fco/2fc^)]- 



:^Ao- 



[ln(PrV2)] 



(49) 



Note that the particular value of ko is unimportant here, because ln(A:o/2A;^) = ln(fcy/A;^) + ln(fco/2A;i,), the 
second term being subdominant. The above value of A specifies a particular solution of the form (38) which 
exactly vanishes at fc = fco. In §2.3 and §2.4, we saw however that any excitation present at A; > fco at a 
given time must necessarily end up spreading over to the integral scales. The question remains whether, 
once this happens, the entire spectrum would eventually be affected. It is, however, encouraging to notice 
that A = Ao represents a nontrivial limit when the infrared cutoff ko is taken to zero. 
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3. TWO-POINT CORRELATION FUNCTIONS OF THE MAGNETIC FIELD IN THE 

CONFIGURATION SPACE 

We would now like to put the solution of the kinematic- dynanio problem in a somewhat different 
perspective and develop an adequate theory valid in the long-time limit. In order to make any progress, we 
necessarily have to renounce the small-scale approximation. While it has furnished us with a fair amount of 
practical understanding of the initial behavior of the magnetic fluctuations in the kinematic regime, it does 
not contain sufScient information to uniquely fix the growth rates and the spectra eventually attained by 
the magnetic energy spectrum. We must therefore attack the kinematic dynamo problem in the form that 
remains valid at all scales. This can be more conveniently accomplished in the configuration space, for the 
x-space analog (inverse Fourier transform) of the mode-coupling equation (11) is local and can be reduced to 
the quantum mechanics of a particle with variable mass in a one-dimensional potential well. This approach, 
first proposed by Kazantsev (1968), leads to a clearer mathematical and physical picture of the properties 
of the magnetic-field correlation functions in the long-time limit. 

Kazantsev's quantum mechanics was studied in the case of 3D incompressible fiow by several authors 
in 1980s (Ruzmaikin and Sokoloff 1981; Novikov, Ruzmaikin, and Sokoloff 1983; Artamonova and Sokoloff 
1986; Klceorin, Ruzmaikin, and Sokoloff 1986; Maslova and Ruzmaikin 1987), mostly for particular forms of 
the velocity correlation function. Novikov, Ruzmaikin, and Sokoloff (1983) also set forth several important 
ideas with regard to the 2D case. In this Section, we develop the configuration-space two-point kinematic 
dynamo theory for a general, d-dimensional and arbitrarily compressible advecting flow. In § 3.1, we derive an 
evolution equation for the correlation function of the magnetic field. This equation is the configuration-space 
analog of the mode-coupling equation of the previous Section. We then proceed to formulate the Kazantsev 
quantum- mechanical form of this equation. In § 3.2, the three principal asymptotic regimes of this model 
are considered: the asymptotic solutions are found in the subresistive, inertial, and integral scale ranges. 
In § 3.3, these solutions are matched in a systematic way and the dynamo growth rates are determined for a 
particular choice of small- and large-scale regularization of the Kazantsev quantum mechanics. The ensuing 
results have universal applicability if the specific form of the regularization is unimportant. The possible 
relevance of the large-scale structure of the velocity correlations is further discussed in § 3.4. 



3.1. The Kgizantsev Equation 

We start with the induction equation in its standard configuration-space form (1). The usual assump- 
tions of homogeneity, isotropy, and mirror invariance, as well as the Gaussian white-noise nature of the 
velocity field are made [Eq. (2)]. The velocity correlation tensor in the x space can be represented as follows: 

.-(y) = ^L.iy)S- - (...(.) - ^) - ^) , (50) 

where KLL{y) and K^Niv) are called the longitudinal and the transverse correlation functions. The ten- 
sor K*-' (x) defined above is the Fourier transform of the k-space correlation tensor k^^ (k) whith which we 
operated in the preceding Section [see Eq. (5)]. The exact relationship between the functions KLL{y), K,NN{y) 
and the k-space correlation functions K{k), k{k) is worked out in detail in Appendix A.l. 

We now proceed to derive an evolution equation for the correlation tensor H^^{t,y) of the magnetic 
field [see definition (3)] by directly averaging the dynamic equation for the tensor B^{t,x)B^ {t,x') (here, as 
always, x — x' = y) and splitting the arising triple averages according to the Furutsu-Novikov formula (see 
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Appendix B). Namely, 

_d 
dy 

where the four terms on the right-hand side are the triple averages: 



dtW'iy) - 2r]AW^{y) = ^ C*^ (y) - C'^'^y) - C^^\-y) + C^--"(-y) , (51) 



(7«^ (y) = C"=^^-(xi - X2) = (w'=(xi)B^(xi)B^(x2)) 

= I [(«"(y) - «'='(0))i?,7(y) - 4^(y)i/''(y) + ^'H^'iy)] ■ (52) 

The details of the derivation of the above expression are standard and straightforward and have therefore 
been omitted. The reader is reminded that lower indices occurring after a comma designate derivatives with 
respect to y. Upon collecting all terms in Eq. (51), we arrive at the following equation for the correlation 
tensor H^^{t,y) (of. Vainshtcin 1980, 1982b; Molchanov, Ruzmaikin, and Sokoloff 1983; Kichatinov 1985; 
Vainshtein and Kichatinov 1986): 

-2K%Hy - K%iW^ + ,^^iH^i + K%H^\ (53) 

Since the magnetic field is isotropic and solenoidal (iJ^ = 0), all the necessary information about its 
two-point correlation properties is contained in just one scalar function of y. The most convenient (and 
customary) choice is the longitudinal correlation function 

HLL{t,y) = ^H^\t,y), (54) 

in terms of which the correlation tensor W^{t, y) can be expressed as follows (see Appendix A. 2): 

H'i {t, y) = HLLit, y)6'^ + ^ yH'^Lit, v) (s'^ - ^) • (55) 

Here and below primes denote derivatives with respect to y. Multiplying Eq. (53) by y^y^ /y'^ and making 
use of the formulas (54) and (55), and after more tedious algebra, we find that the longitudinal correlation 
function of the magnetic field satisfies the following evolution equation: 

+ ^ {K'{y) + Q'{y) + ^ [K{y) - Q{y)\ ) H^l, (56) 



where we have introduced two "rcnormalizod diffusivit 



ICS 



K{y) = 2r) + KLL{Q)-nLL{y). (57) 

= 2, + : I'-^y^ = 2, + .,,(0) - (58) 

Note that, in the incompressible case, we should have (see Appendix A. 2) 

Q{y) = K{y) + y^. (59) 
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Due to the realizability constraints, it is assured that K{y) > 0, Q{y) > (sex; Appendix A. 3). Eq. (56) is the 
configuration-space analog of the mode-couphng equation (11). While we have chosen to derive it by directly 
averaging the induction equation, it could also have been done by Fourier-transforming the mode-coupling 
equation (11). We should like to remark, however, that Eq. (56) is more convenient than the mode-coupling 
equation (11) for the study of general velocity correlators, such as those whose small-scale expansion is 
nonanalytic (as in the case of Pr < 1). 

A further modification of Eq. (56) brings us to the quantum-mechanical analogy that was mentioned in 
the Introduction. Namely, we introduce the "-^-function" 



y'^-'VmHLL{t,y)e.p J'dy'^L^ (60) 



and find that it satisfies the one-dimensional Schrodingcr equation, 

atV = ^;^V'"-^^(y)V', y>o, (61) 

which describes the quantum mechanics of a particle with variable mass^ miy) = [2K{y)]~^ in a one- 
dimensional (or radial) potential 

Eq. (61) [with potential (62)] was first derived by Boldyrev and Schekochihin (1998). In the three-dimensional 
case and for incompressible velocity field, it reduces to the equation derived by Kazantscv (1968) (sec also 
Kichatinov 1985; Vainshtein and Kichatinov 1986). In what follows, the general equation (61) will also 
be referred to as the Kazantsev equation. The expression for the potential (62) in the incompressible case 
becomes [see formula (59)] 



(63) 



It is this incompressible form that has mostly been considered in the literature (Ruzmaikin and Sokoloff 

1981; Novikov, Ruzmaikin, and Sokoloff 1983: Artamonova and Sokoloff 1986; Kleeorin, Ruzmaikin, and 
Sokoloff 1986; Maslova and Ruzmaikin 1987). Note that, in 2D incompressible case, it is immediately clear 
from the formula (63) that the potential is repelling and hence no bound states exist. 



3.2. The Eigenvalue Problem Associated with the Kinematic Dynamo: Asymptotic 

Solutions 

We now consider the eigenvalue problem associated with the Kazantsev equation: 

-\7^lj = ~^^^r + V{y)il;, y>0, (64) 

where, for the sake of theoretical uniformity, we measure the eigenvalues (dynamo growth rates) in terms 
of 7, the magnetic energy growth rate (14) that figured in Sec. 2. [In terms of the configuration-space 



^The mass m{y) is always positive due to the realizability constraint (A38) (see Appendix A.3). 
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quantities. 7 is proportional to the square of the velocity shear: 7 = |k^^(0)|/2^, where A is defined by the 
formula (22). See Appendix A. 4, formula (A51).] The objective is to find the spectrum of the Schrodinger 
operator in the right-hand side of Eq. (64). Clearly, positive energies correspond to damped modes (A < 0), 
while negative energies represent growth (the dynamo effect, A > 0). Therefore, the ground state of Eq. (64), 
if it exists, would give the desired long-time asymptotic of the magnetic- field correlation function. 

At small scales, the velocity correlation tensor (50) has a Taylor expansion, 

(y) = - 1 {y^6'^ + 2a fy^) + • • • , (65) 

and the renormalized diffusivities K{y) and Q{y) therefore are, to second order (see Appendix A. 4), 

K{y) ~ 2r,+ i(l + 2a)«2y', (66) 
Q{y) ~ 2T]+^K2y\ (67) 

where a is the compressibility parameter. These are valid for y <C («;o/'«2)^^^ ^ k~^, where kq = kll{0)- 
Note that (1 -|- 2a)K2 = \k'Ii^{0)\ = 2jA [see formula (A51)]. 

When y (ko/k;2)^^^, both K{y) and Q{y) tend to a constant value: 

K{y) ^ Q{y) ~ 2?? + kq ^ kq. (68) 



Substituting these small- and large-scale asymptotics of K{y) and Q{y) into the expressions for the 
mass of our particle and for the potential it lives in [see (62)], we may gain a qualitative idea of what m{y) 
and V{y) look like. [The typical forms of the mass and potential (in 3D) are sketched in Fig. 7 and Fig. 8.] 
Three regions with diflferent asymptotics are clearly pronounced: the subresistive scales, the integral scales, 
and the intermediate inertial range. Let us analyze the situation in each of them in turn. 



The subresistive scales. These are, of course, the scales satisfying y <C (??/k2)^^^ ~ kn^- There, the 
mass of the particle is constant and the potential is a repelling inverse square: 

m(y) ~ ^, (69) 
The asymptotic solutions of the eigenvalue problem (64) in this approximation are 

Hy) = y^'^'^/', (71) 

where it is clearly the "-|-" -branch that must be picked. Indeed, it satisfies just the right boundary condition 
that follows from the definition of the "wave function" (60) and the requirement that HLL{t, y) have a positive 
finite value at y = (it is the total magnetic energy in the system). Note that the boundary condition could 

not have been satisfied had there not been a nonzero magnetic dijfusivity in the problem. 



The solution in this region does not depend on the value of A, in a precise analogy to what we saw in 
the wave- number space [see the asymptotic (45)]. 
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The integral scales. Considering the opposite end of the available range of scales, for y [kq/kq)^/'^ ~ 
k~^, we again find a constant (now "renormalized" ) mass and a repelling inverse-square potential: 

This time, however, it is the potential that contributes subdominantly to the solution, while the eigenvalue A 
plays an important role (unless A = 0). The solutions that do not grow at infinity are 

^{y) = exp(^-^y^, A>0, (74) 
m = y-^'-'^/', A = 0, (75) 



i;{y) = sin|^^^y + ^j, A < 0. (76) 

We see from the form of the solution, as well as from the fact that the potential approaches zero from above 
as y +00, that no bound states can exist with negative A. In fact, neither can they exist with A = 0, for it 
leads to a power-like decay of the wave function at best. 



The inertial range. These are the scales that are equally removed both from where resistive-diffusive 
effects operate and from the integral scales where the details of the velocity-field correlation properties start 
to be felt. Namely, {r]/K2y^^ <C y <C («;o/'«2)^^^. In this regime, we obtain 



V{y) ~ -Vo, where Vq = ^ 



2(1 + 2a) 



(78) 



We have clearly found ourselves inside the potential well, which is locally flat. An obvious estimate for the 
eigenvalue immediately follows: A < Vb/7 = Amax, which is the same as the upper bound obtained in §2.5 
[sec the inequality (43)]. We find that, in d = 3, the potential well always extends below zero (so Vq > 0), 
while in d = 2, < for a < -1/4. 

Eq. (64) in this regime takes the form 

V'" + ^^^i^V'-0, (79) 

where A is defined by the formula (22) and it can be easily verified that (1 + 2a)K2 = 2'jA. The effective 
potential in Eq. (79) is of the inverse-square kind, as advertised in the Introduction (§ 1.3). This equation is, 
of course, the configuration-space version of Eq. (40) in the inertial-scale limit. Its solutions are, expectedly, 

V'(y) = Ci2/l/2+V(A-Ao)M + c2j/l/2-V(A-Ao)M^ ^gg) 

where Aq = Amax — A/4: is defined as it was in § 2.5 [see formula (28)]. 
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3.3. Asymptotic Matching: The Dynamo Growth Rates in the Long-Time Limit 

In order to complete the asymptotic solution of the problem, one must now proceed to properly match 
the solutions we have obtained in the three asymptotic scale ranges described above. 

Matching the inertial range solution with the subresistive one does not present any difficulty in principle. 

Indeed, wc saw in Sec. 2 that, in the k space, an analytic sohition could be obtained that is valid both in the 
resistive and inertial ranges. The issue of matching the inertial-range and the large-scale solutions is more 
troublesome. The appropriate matching procedure depends on the full functional form of the potential at 
the velocity-field scales and must typically be performed assuming a particular form of the velocity correlator 
(cf. Artamonova and Sokoloff 1986; Kleeorin, Ruzmaikin, and Sokoloff 1986). 

In this section, wc propose a simple way to match the asymptotic solutions obtained in § 3.2 and 
calculate the dynamo growth rates in the long-time limit. Our method is based on the assumption that 
the particular way the Kazantsev particle mass and potential are regularized at small and large scales does 
not play an essential part in determining the dynamo growth rates. This assumption holds for such forms 
of the velocity correlation function that ensure that the potential well be located at small (inertial) scales 
and that the potential (62) have no relevant structure at the integral scales. In other words, our model 
quantum-mechanical particle must be "localized" at small scales in order for the structure of the velocity 
correlations to be unimportant. If this is satisfied, the dynamo growth rates in the long-time limit can be 
calculated in a universal way. Let us explain how this can be accomplished. 

Although, as we have noted, matching the subresistive and the inertial-range solutions is no problem in 
principle, it does require a considerable amount of algebra (Kleeorin, Ruzmaikin, and Sokoloff 1986). The 
formal difficulty lies in the fact that tp{y) and tp'{y) at y ~ (?7/«;2)^^^ are of different orders in rj. The 
matching must therefore be carried out in two successive orders of the asymptotic expansion. This difficulty 
can be circumvented in the following way, suggested by Landau's discussion of the quantum mechanics in 
the inverse-square potential (Landau and Lifschitz 1977). Assuming that the particular form of the small- 
scale regularization is not essential for choosing the right solution in the inertial range, let us pick a small 
Vv ^ iv/f^2y^'^, extend the inertial-range forms of the mass and potential [formulas (77) and (78)] down to 
this point (to all y > yr,), and replace the actual values oim{y) and V{y) at y < by constants equal to their 
values a.t y = yrj (thus taking "a shortcut" to the "wall" at y = 0). Let us then solve the resulting eigenvalue 
problem in both regions and match the solutions at y = y^. We will then carry out a similar procedure to 
match the inertial-range and the large-scale solutions at some yo ^ (ko/k2)^^^ (see below). Eventually, we 
must take the limit yr, ^ and yo ^ oo and determine the eigenvalues A from the requirement that the 
matching relations at yr, and yo be consistent with each other in this limit. 

Let us empahsize that our matching procedure corresponds to a particular reasonable choice of the 
small- and large-scale regularization of the Kazantsev quantum mechanics. Its universal applicapility hinges 
on the assumption that the specific form of such regularization is inessential. 



The proper solution in the interval < y < y,, is ipiu) — sin(7,,y/r/,,), where 7^ = 1/ (A„iax — ^)/^ (we 
assume, of course, that A < Amax)- The solution at y > y,, is given by the formula (80). Matching tp{yn) and 
V''(yr,), we get, for A > Ao, 

c^^s^-jr^ Ao<A<A_. (81) 

ci 'Yn cot Jrj- S2 ' 



Here we have denoted si,2 = 1/2 ± ^(A — Ao)/A. For A < Aq, it is convenient to write the inertial-range 
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solution (80) as 



ip{y) = cy^^"^ sin 



5X = Xq- X>0. 



The matching determines the value of the phase: 

_ 7^ cot 7^ - 1/2 



cotx = 



where 7-7= 2^/1 + 



(82) 



(83) 



Let us now employ the same device to match the inertial-range and the large-scale solutions. Namely, 
we now pick some yo ~ (^0/^2)^^^, formally extend the validity of the incrtial-range asymptotics of m,{y) 
and V{y) to all y < yo, while at y > yo, let m{y) = 1/2ko and V{y) = 0. The solution in the latter region 
is the large-scale solution (74) or (76) (depending on the sign of A). For convenience, let us take yo = 
{Ko/'yAy/'^ , whence, in the formulas (74) and (76), wc have, by definition, X^/kq = X/Ayfj = (7o/yo)^j 
where 70 = yjAjTA. In order to match tp{yo) and ip'{yo), we take the asymptotic solutions appropriate for 
various values of A: for A > Aq, the formulas (80), (74), and (76) imply 



Cl 
C2 

£1 

C2 



S2 



„«2-Sl 



Si +70 

S2 - 70 cot (70 + ip) 



A > Ac & A > 0, 



7ocot(7o + ip) - si 
for A < Ac, we use the formula (82) instead of (80) and get 



Ao<A<0; 



cot 



cot 



'^4 ln(^ 
'4 In ( ^ 



X 



X 



70 + 1/2 

^/5X/A ' 

70 cot (70 + ip) -1/2 



^/5X/A 



< A < Ao, 

A < Ao & A < 0, 



(84) 
(85) 

(86) 
(87) 



where ^A = Aq — A. Note that the relations (84) and (86) smoothly extend to the case of A = (where 
7o = 0), which does not have to be treated separately. 

We are now in a position to find the dynamo growth rates. In the matching conditions derived above, 
we take the limit y^/yo and determine the allowed values of A from the consistency requirements. 



The case A > Aq & A > (growth). For these values of A, the relations (81) and (84) that ensure proper 
matching of the inertial-range solution at dissipative and integral scales respectively are incompatible in the 
limit yn/yo — * unless 

(70 + si ) (7,, cot 7,, - S2) = 0. (88) 

This is a transcendental equation on A, which is readily seen to have no solutions. Thus, no positive growth 
rates that are larger than A07 are allowed. 

Note that, in 2D, while we have Ao < 0, a rudimentary potential well exists for some values of the 
compressibility parameter: in §3.2, we already pointed out that Amax > for a > —1/4. We see now that 
this well has turned out to be too shallow to sustain a bound state with A > > Aq. 
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The case < A < Aq (growth). In this case, we have to make sure that the relations (83) and (86) are 
consistent with each other in the limit yo/z/?? ~ Pr^^^ ^ cxd. If (5A = A — Aq remains finite under this limit, 
the left-hand side of Eq. (86) has no definite limit, while its right-hand side is equal to a finite constant. 
The matching is therefore impossible. On the other hand, if we assume that SX — > -|-0 as yo/Vr) — * oo, the 
right-hand sides of both relations (86) and (83) tend to infinity, whence x = ""Tn, m € Z, and 

^A= ^~ — -^<1, n = l,2,... (89) 

Clearly, it docs not matter here what the particular definitions of j/o and are: miiltiplying any given such 
definitions by finite constants leads to a finite, and hence subdominant, correction in the denominator of the 
expression (89). 

The expression (89) reproduces the result that was obtained in § 2.5 [see formula (49)] by requiring the 
magnetic-energy spectrum to vanish at some particular ("smallish") wave number. Within our model, we 
have now demonstrated that this result gives a valid long-time asymptotic for the kinematic dynamo problem. 
The key to understanding this coincidence lies in the fact that the magnetic-field correlation function decays 
exponentially at large scales [asymptotic solution (74)]. Due to the solenoidality of the magnetic field, such 
exponential tail implies that -ff(fc) must vanish at fc = [see Appendix A. 2, formula (A27)]. Imposing this 
as a boundary condition on the small-scale spectrum (44) necessarily implies (49) or (89). 

The solution we have thus obtained is relevant in three dimensions, where Aq > 0. Note that formula (89) 
can be used to estimate the critical value of the Prandtl number that has to be exceeded in order for the 
dynamo effect to occur (i.e., for the growth rate to be positive, A = Aq — (5A > 0): Pre exp(27ri/j4/Ao). 
For the incompressible case, Pre ^ 26; for the irrotational case. Pre ^ 17000. It must be clear, of course, 
that such critical values are not to be taken as anything more than rough estimates, for they have been 
deduced from the formula (89), which is strictly valid only when ^A ^ 1. Numerical studies of the excitation 
threshold for the Kazantsev equation give (in the incompressible case) values of the critical Prandtl number 
(or critical magnetic Reynolds number, which, for a one-scale velocity field, is the same thing) about Pre ^ 53 
(Novikov, Ruzmaikin, and Sokoloff 1983; Maslova and Ruzmaikin 1987). 



The case Aq < A < (damping). Here we must reconcile the relations (85) and (81) in the limit yr^/yo 
0. In view of the incompatibility of Eq. (88), this can only be done if 

70 cot(7o -I- <^) - si = 0. (90) 

This fixes the value of the phase for any applicable A. No further constraint ensues for A lying in the 
interval Aq < A < 0. Provided this interval is nonempty (as is the case in two dimensions, unless a = 0), a 
continuous spectrum of eigenvalues exists and fills the entire interval. Since the upper bound of the spectrum 
is zero, a power-like decay of the magnetic fiuctuations should be expected. 



The case A < Aq & A < (damping). In this case, just like in the case of < A < Aq, assuming that 5X 

remains finite (or becomes large) a.s yo/yrj ^ oo leads to incompatibility between the right- and left-hand 
sides of Eq. (87): on the right-hand side, there must be a finite constant, while the left-hand side has no 
definite limit. If Aq > (3D), this is the only possibility. Thus, if Xq is positive, negative values of X are 
ruled out. 

If Ao < (2D), it is possible to have ^A +0 as yo/yri — * oo- Again, the argument that led to the 
solution (89) remains valid. However, this solution is not relevant because we have already discovered a 
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continuous spectrum between Aq and 0. It is the resulting power- like decay that will determine the long-time 
behavior of the magnetic fluctuations. The case of Aq = 0, which occurs in 2D for a = (the "quasiequi- 
librium" case we noted in §2.3), appears to be exceptional in this context, for then the interval (Ao,0) 
is empty and the only applicable solution is, in fact, (89) (where now A < Aq = 0). Namely, we have 
A ~ — ^TT^/ [ln(Pr^/^)] —0, which should also result in a power-like temporal decay law in the limit of 
vanishing resistivity {rj +0). 

Remark on the case A = Aq. When A = Aq, we have Si = S2 = 1/2 and the solution (80) appears 
to be degenerate. This degeneracy is spurious. In fact, the power-like solutions (80) are asymptotics of 
the modified Bessel functions / 77- — , , , . and K n- — , , , , . This can be seen if the Kazantsev eigenvalue 

problem (64) is solved in the incrtial range with next-order correction terms in the Taylor expansions of K{y) 
and Q{y) included [see Eq. (C8) in Appendix C below]. When A = Ao, we should therefore use, instead of 
the asymptotic solution (80), 

My)=y'^Hci+C2liiy). (91) 

The matching at t/^ and j/o can be carried out in quite the same fashion as before. The consistency conditions 
as ln(yo/y»)) 00 are given by Eq. (88) if Aq > (still incompatible) and Eq. (90) if Aq < (can always be 
satisfied). 



3.4. Discussion 

Thus, assuming that the particular form of the large-scale rcgularization of the Kazantsev quantum 
mechanics is unimportant, or, in other words, that the Kazantsev quantum-mechanical particle is "localized" 
at small (inertial) scales, we have been able to make progress in our investigation of the long-time asymptotic 
behavior of the passive magnetic fluctuations. In 2D, a continuous spectrum of negative A's exists, which 
indicates a power-like temporal decay. Its exact nature can be discovered only by solving the time-dependent 
dynamo equation (56). In 3D, a universal dynamo growth rate has been found, which asymptotically 
approaches A07 from below, as Pr — > 00 (the convergence is square-logarithmic). 

The exact quantitative criterion that would clearly define the class of velocity correlators within which 
such universality with respect to the large-scale rcgularization exists, is unknown. Due to the complicated 
structure of the potential (62), it appears to be very hard to determine. On the other hand, no examples 
of realizable^ velocity correlators for which the dynamo growth rate is essentially determined by the large- 
scale properties of the velocity field are available. As we already mentioned in the Introduction, numerical 
solution of the Kazantsev eigenvalue problem (64) for various plausible forms of the velocity correlator 
(see Novikov, Ruzmaikin, and Sokoloff 1983; Artamonova and Sokoloff 1986; Maslova and Ruzmaikin 1987, 
and Appendix C), including the one corresponding to the Kolmogorov turbulent spectrum (Kulsrud and 
Anderson 1992), did not reveal any nonuniversal behavior of the second-order statistics of the magnetic field. 
Thus, the validity of the universal results presented in this section appears quite robust. 

A tentative explanation for such robustness could be given along the following lines. The mass mly) 
of the Kazantsev particle at integral scales tends to the constant finite value m(oo) ~ 1/2ko. At small 
(inertial) scales, however, it sharply increases and, at the resistive scales, reaches the asymptotically infinite 



''I.e., those whose Fourier transforms are correlation functions (see Appendix A.l). 
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value m(0) = l/4?7 (for -q +0). Thus, the particle is impelled to slide toward the small scales not just by 
the potential, but also by its own variable mass: it is "heavier at small scales." This leads to the large-scale 
effects being marginalized even if the value of the potential in the inertial range, V{y) ~ —Vq, is not, in 
fact, its global (or even local) minimum. Note that, while the latter situation has not been encountered in 
the previous studies of the incompressible case, it turns out to be quite common when the advecting flow 
possesses some degree of compressibility (see Appendix C). 

It remains to be seen if examples of velocity correlators can be constructed for which the second-order 
statistics depart from the universal solutions obtained above. The issue of whether higher-order statistics 
may be more vulnerable in this respect is also open to investigation. At the current level of understanding of 
the two-point kinematic dynamo, we may conclude that, while the nonuniversality of the long-time solutions 
due to large-scale efi'ects remains no less of a possibility than it was when first hinted at by Kazantsev (1968), 
the results obtained assuming small-scale universality appear quite reliable for practical purposes. 

4. Where Do We Go From Here? 

The main physical conclusion of the kinematic theory developed above is that, at the moment when 
the nonlinear effects become important, the magnetic energy is concentrated at the smallest scales that it is 
able to reach at the exponential rate ~ 7 during the kinematic stage of the dynamo. An ascending power 
spectrum k^/^ extends through the magnetic inertial range up to the wavenumber where the peak of the 
spectrum is located at the time the nonlinearity sets in. 

It is important to understand that the second-order (spectral) theory developed here docs not contain all 
of the statistical information necessary to determine the condition for the onset of the nonlinearity. Indeed, 
the nonlinear back-reaction effect is controlled by the Lorentz tension force oc (B • V)B in the momentum 
equation of the fluid (in the incompressible case). The second-order statistics of this quantity are fourth order 
in the magnetic field. Importantly, the statistics of (B • V)B arc determined not just by the magnitude of 
the magnetic field and by its overall characteristic scale, but also by the structure of the field, namely, by the 
characteristic scale of its variation along itself. This issue was thoroughly investigated by Ott and coworkers 
(see review by Ott 1998), Cattaneo (1997, 1999), and by Schekochihin et al. (2002a) (see also references 
therein). The result is that the parallel scale does not decay with the overall characteristic scale, but stays 
at values comparable to the velocity scale k~^. The small-scale field is arranged in a pattern of randomly 
oriented long (of characteristic size ~ fc~^) folds, whithin which the field undergoes rapid transverse spatial 
oscillations of its direction. The implication is that the condition for the onset of the nonlinearity is the 
approximate equalization of the magnetic- field energy and the energy of the smallest turbulent eddies. 

Once this has occurred, a nonlinear theory is needed to describe futhcr evolution of the field. No such 
theory exists as yet. However, there does exist a fair amount of numerical evidence, which, along with the 
heuristic physical insight, can be used as guidance in mapping out the tentative strategies of picking up 
where the kinematic theory left off and building on its results. Resolving three main issues appears to be 
the first order of business at this point. 

Nonlinecir saturation at subviscous scales. The onset of the nonlinear regime can occur under two 
distinct sets of circumstances, depending on the magnitude of the initial seed field and on the magnetic 
Prandtl number. As we saw in the preceding sections, the growth of the field strength and the refinement 
of its scale proceed in parallel. With the seed field assumed to be concentrated at velocity scales ^ the 
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time it takes for the dynamo to switch from the diffusion-free to the resistive regime is of order ~ 7"^ logPr. 
Given a very weak seed field, the dynamo will still be kinematic at the time of this transition. The input 
state for the nonlinear theory will then be that with the small-scale magnetic-energy spectrum extending 
through the entire subviscous scale range and peaked at the resistive scale ^ (see § 2.4, § 2.5 and Sec. 3). 
Alternatively, it is conceivable that the initial field should be sufficiently strong for the nonlinear regime 
to commence before the resistive scale is reached (Kulsrud 2000) (see §2.3). One wonders then whether 
the nonlinear back reaction is capable of arresting all further refinement of the field scale. Although the 
answer to this question is in the negative, it appears that the spreading of the magnetic excitation into 
the remainder of the subviscous range available to it proceeds at a much slower pace than in the kinematic 
regime: most likely, at the resistive time scale (Kinney et al. 2000a,b; Schekochihin et al. 2002b). It may be 
possible to study the evolution of the magnetic spectrum through this initial nonlinear stage in a theoretical 
framework that is substantially simpler than the full nonlinear MHD. Namely, the action of the fluid motions 
on the small-scale magnetic fields can be modelled by balancing the viscous and magnetic stresses, rather 
than by solving the full Navier-Stokes equation with Lorentz back reaction (the so-called Stokes model). 
An "intermediate" small-scale nonlinear theory emerges from such an approach that is a direct extension of 
the kinematic theory elaborated above (Kinney et al. 2000a; Schekochihin et al. 2002b); an exactly solvable 
model of such a regime (without diffusion) was constructed by Boldyrev (2001). 

Plasma damping mechanisms. A parallel question is how well the Spitzer resistive-diffusive small- 
scale regularization that is a part of the standard MHD paradigm describes what really happens at subvis- 
cous scales in the high-Prandtl-mimber astrophysical plasmas. It must be realized that, while the resistive 
scales fc~^ that obtain by substituting the ISM or protogalactic- plasma parameters into the Spitzer formula 
for the magnetic diffusivity r] may be 7 to 11 decades below the viscous scale k~^, the mean free path of the 
particles that make up these plasmas is substantially larger (the simplest estimate gives ^mfp '"^ ^ kj^ ^ , 
i.e., only about 10 times shorter than the viscous scale). The MHD approximation is not strictly valid 
at scales below the mean free path, and one must allow that a number of plasma damping mechanisms 
may supercede the Ohmic resistive dissipation (Vainshtein 1976; Kulsrud and Anderson 1992; Kulsrud et al. 
1997b). These include, most importantly, the ambipolar damping that is present due to the partial ionization 
of the ISM (Kulsrud and Anderson 1992; Chandran 1998; Subramanian 1997, 1998, 1999; Brandenburg and 
Subramanian 2000) as well as the kinetic effects due to the Braginskii (1965a) tensor viscosity and magnetic 
diffusivity (Malyshkin 2001) and to Landau damping of magnetic fluctuations (Kulsrud et al. 1997b). We 
note that, while all of these effects are nonlinear in the magnetic field, the kinetic effects depend solely on the 
direction, not the magnitude, of the field, which only needs to be sufficient to magnetize the ions — a condi- 
tion that is satisfied in the astrophysical plasmas in question already for very weak fields, which are otherwise 
passive. Treatment of these effects again appears possible by means of certain amendments to the kinematic 
theory developed here (cf. Kulsrud and Anderson 1992; Subramanian 1997, 1998, 1999; Schekochihin and 
Kulsrud 1998; Malyshkin 2001). Note that, in view of the possibility of other-than-resistive small-scale reg- 
ularization, the study of the onset of back reaction in the diffusion-free dynamo, as discussed in the previous 
paragraph, acquires additional importance. 

Growth of large-scale magnetic fields. Finally, the most important — and difficult to resolve — issue 
is that of the feasibility of some form of nonlinear transfer of the magnetic energy accumulated at small 
scales toward large scales: the inverse cascade. Two distinct directions of inquiry must be clearly identified 
here. 
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First, as wc mentioned in § 1.1, it has long been known that the presence of a helical component in the 
fluid turbulence leads to exponential growth of the large-scale magnetic field in the kinematic regime {the 
a effect). In the high-Prandtl- number astrophysical plasmas, this growth is much slower than that of the 
small-scale fluctuations. Once the latter have entered the nonlinear regime, they are expected to substantially 
modify the a effect, almost certainly suppressing it to some extent. The magnitude of this suppression is 
a subject of an ongoing debate, various estimates ranging from extreme quenching with suppression factors 
as large as Rm to moderate order-one reductions (see, e.g, discussion in a recent series of papers by Field, 
Blackman, and Chou 1999; Blackman and Field 1999, 2000; Chou 2001a; Field and Blackman 2001, and 
references therein). While none of the existing theories of the nonlinear a effect is universally agreed on, it 
appears safe to expect that (i) helicity does lead to some growth of the large-scale magnetic field; (ii) this 
growth is slowed down in a nonlinear fashion by the accumulating small-scale magnetic fluctuations. At large 
Prandtl numbers, the suppression effect may be very large, making impossible any growth of the large-scale 
fields at time scales faster than the resistive. The kinematic a-effcct growth rate, which is already about 
10000 times smaller than that of the small-scale dynamo, is thus further reduced by a factor possibly as large 
as Km- In the context of the ISM's nearly perfect conductivity, such a dynamo would be quite prohibitively 
slow and could not explain the observed galactic magnetic field. [For some numerical evidence regarding 
these processes, we address the reader to papers by Meneguzzi, Frisch, and Pouquet (1981); Cattaneo and 
Hughes (1996); Brandenburg (2001); Maron and Blackman (2001). However, none of these simulations has 
enough resolution to be conclusive with regard to the high-Prandtl- number plasmas.] 

The other avenue of investigation, which we announced in § 1.1, is to look for the possibility of nonhelical 
inverse cascade leading to eventual equipartition of the magnetic and kinetic energies at all scales including 
the large. Indeed, such an eventual equipartition appears eminently plausible at the flrst glance, since it 
corresponds to the steady-state spectra that follow from Kolmogorov-style theories based on the view of the 
MHD turbulence as resulting from nonlinear interactions of Alfven waves (Iroshnikov 1964; Kraichnan 1965; 
Goldreich and Sridhar 1995, 1997, — in this context, the outcome of the ongoing controversy over the exact 
form of the saturated equipartition spectrum is immaterial). However, several recent numerical simulations 
with Pr > 1 and no externally imposed dc background fleld have cast serious doubt on the possibility of the 
saturated equipartition state (Kinney et al. 2000a, b; Chou 2001b; Maron, Cowley, and McWilliams 2002; 
Schekochihin et al. 2002b). The questions of what the saturated magnetic spectrum is and whether the 
magnetic energy at the large scales can be shown to achieve values sufficient to give rise to an adequate 
theory of the observed galactic field, remain unresolved. 

We are grateful to J. A. Krommes who read an earlier manuscript of this work. Both the substance and 
the style of the presentation have benefited from his suggestions. We would also like to thank E. Blackman, 
S. Cowley, J. Maron, V. Pariev, V. Rytchkov, D. Uzdensky, and the anonymous referee for useful comments. 
This work was supported in part by the U. S. Department of Energy under Contract No. DE-AC02-76-CHO- 
3073. 



A. CORRELATION FUNCTIONS OF ISOTROPIC RANDOM FIELDS AND OTHER 

MATTERS 

In this Appendix, we provide a set of useful facts and formulas that are frequently needed in turbulence 
calculations such as those of Section 2 and Section 3. We will work with a random field u% which we will 
often refer to as the velocity field. However, all the results of this Appendix also hold for the magnetic 
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field provided solenoidality constraint is properly implemented. We find it most convenient to work in 
d dimensions. For d = 3, many of the results below can be found in the lucid and comprehensive chapter on 
the mathematical description of turbulence by Monin and Yaglom (1975). 

This Appendix is organized as follows. In § A.l, we provide the transformation formulas that relate the 
correlation functions of isotropic fields in configuration and wave- number (Fourier) spaces. These impose 
certain constraints on the classes of fimctions in the x space that are proper correlation functions. In § A.2, 
the additional constraints arising from solcnoidal and potential nature of the fields are explained. In § A. 3, 
the so-called realizability conditions are discussed that determine some essential features of the theories 
dealing with isotropic random fields. In §A.4, we give some relations between the small-scale-expansion 
coefficients of the two-point correlation functions in the x space and their counterparts in the k space. 
These relations provide a bridge between the formalisms of Sec. 2 and Sec. 3 and are also needed to put the 
results of this paper and the one-point statistical results published earlier (Boldyrev and Schekochihin 2000; 
Schekochihin and Kulsrud 2001; Boldyrev 2000, 2001; Schekochihin et al. 2002a) in a unified theoretical 
framework. Finally, in §A.5, we list the angle-integration identities that were used in the derivations of 
Sec. 2. 



A.l. Isotropic Fourier Treinsforms 

We will employ the usual Fourier- transform conventions: 

u\k) = I d'^xe-'^-''u'{-K.) (direct), (Al) 



In this Appendix, we are not concerned with time dependence of our fields and concentrate purely on their 
properties arising from the assumptions of spatial homogeneity, isotropy, and mirror invariance. Under this 
set of assumptions, we may define the correlation functions of the field as follows: in the configuration 
space, 

(tt^(x)u^(x')) = K*^(x-x'), (A3) 



I . KNN{y) 



- ^ ) , (A4) 



where KLhiy) and KNN{y) are called the longitudinal and the transverse correlation functions for obvious 
reasons. In the Fourier space, 

(u'(k)MJ(k')) = {2'K)'^5{k + k')K'^{k), (A5) 

K'^{k) = K{k)S'^ +K{k)'^. (A6) 

The correlation tensors (y) and At*-' (k) defined above are Fourier transforms of one another: 

'^'^■(y)= / ^e^'^-^n'^ik). (A7) 
Let us substitute (A6) into the above formula and find that 

K»i(y) = 6'^ f e^^-y^ik) - / e'^^ y'^. (A8) 
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The d-dimensional Fourier transforms of the isotropic (radial) functions can be found with the aid of the 
following handy formula (see, e.g., Stein 1993): 



/ 



dOd e 



ik-y 



(A9) 



(/jy)(d-2)/2 

where d^a is the d-dimensional angular differential, and Sd = J d^a = /r{d/2) is the area of the unit 
sphere in d dimensions. We have defined the function 

J{d-2)/2{z) 



^d{z) = r{d/2) ■ 



(AlO) 



(^/2)(d-2)/2' 

which is the kernel of the "isotropic Fourier transform" and has many agreeable properties. Two of them 
will be particularly important for us: ^d{z) is analytic, its Taylor expansion for small z being 



Mz - 0) = 1 - - 



2d 8d{d + 2) 

and ^d{z) satisfies the following differential recursion relation: 



ld_ 

z dz 



Mz) = (-1)"2-" 



^d+2n{z), 



(All) 



(A12) 



where {d/2)n = r(n + d/2)/T{d/2) is the Pochhammer symbol. Note also that ^^{z) = Jo{z) and <^s{z) 
{sva.z)/z. 



The isotropic Fourier transform is then 

d'^A; 



/ 



iky 



K{k) = 



Sd 



dk k''-^^d{ky)K{k). 



(A13) 



{2Try ' ' {2Try Jq 

A similarly defined transform sometimes appears in the literature on special functions under the name of 
Bochner transform. After some straightforward manipulations, we obtain [from (A4) and (A8) by way of 
(A12)] 



KLL{y) = 



Sd 



f 

Jo 



{27ry 

KNNiv) = {d-l) 



dkk'^-'^ 

Sd 
{2TrY 



<^d{ky) [nik) + k{k)] — ^d+2{ky)k{k) 



dk k" 



^d{ky)K{k) + - ^d+2{ky)k{k) 



We may also define the energy (spectrum) function 
I{k) = ^"{k) =dK{k) + k{k), 



i{y) = «"(y) = i^LL{y) + nNN{y) = ^ | dkk''-^^d{ky)i{k). 



(A14) 
(A15) 

(A16) 
(A17) 



The importance of such formulas as (A14), (A15), and (A17), lies in that they impose constraints on 
what functions in the configuration space may in fact be proper correlation functions. Thus, for example, 
since the k-spacc spectrum I{k) must clearly be positive for all k and be properly cut off at small scales, 
formula (A17) stipulates that only an isotropic Fourier transform of such a function can be an x-space 

correlation fimction of w/. 

Note also that the analyticity of the Fourier kernel $d(/cy) [it possesses the Taylor expansion (All) 
around the origin] is the reason for the analyticity of the velocity correlator that passively advected fields 
"feel" at small scales in the Batchelor regime. 
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A.2. Correlation Functions of Solenoidal and Potential Fields 

Two important special cases of the field are incompressible (solenoidal) and irrotational (potential) 
fields. In the Fourier space, we have 

k{k) = -K{k) (solenoidal), (A18) 

K{k) = (potential). (A19) 

In the X space, solenoidality and potentiality impose differential relations between the longitudinal and the 
transverse correlation functions. Namely, 

KNN{y) = yn'^^iy) + {d- l)KLL(y) (solenoidal), (A20) 
K'LLiy) = [y^NNiy) + «ivjv(y)] (potential). (A21) 

[The relation (A20) is known as the von Karman condition.] The equivalence of the above relations in their 
k and x space forms can, of course, be easily demonstrated with the aid of the formulas (A14) and (A15). 

We may summarize the properties of the solenoidal fields as follows 

«^^(y) = >^LL{y)d'' + ^yK'M{s'^-^y (A23) 

where, from (A14) and (A17), 

nLL{y) = 3 77747/ dkk''-^^a+2{ky)m, (A24) 



= ^^A^^(2^)=(^/°°dfc/c^-^$,(Wfc). (A25) 

These formulas are especially important because they must always hold for the correlation functions of the 

magnetic field, which is solenoidal. For example, the formula (A24) imposes a nontrivial (and somewhat 
different from that for the energy) constraint on the class of functions eligible to be longitudinal correlation 
functions of a solenoidal field ■u'(x). 

Note that, if KLL{y) decays at large scales {y oo) faster than l/y"^, then the formula (A25) implies 
an often- used property of the energy function 

POO 

/ dt//-iJ(y)=0. (A26) 
Jo 

This means that I{y) cannot remain positive for all y, and that I{k), which is the inverse Fourier transform 

of /(?/), vanishes at fc = 0. In view of the expansion properties of ^d{ky) [sec (All)], the latter implies that 
the spectrum of a solenoidal field has the following scaling in the limit of small wave numbers: 

fc'^-^7(fc) ~ fc<^+^ as k^O. (A27) 

In 3D, this implies the well-known fc^ infrared scaling of the spectra of solenoidal fields (Monin and Yaglom 
1975). 
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A.3. The Realizability Inequalities 

It is easy to see that the formula (A5), which is a consequence of the homogeneity of space, can be 
recast in the following form: 

K'^ (k) = L-'^ {u\k)u^*{k)) , (A28) 

where L"^ is the volume of space (henceforth set to 1) and the star means complex conjugation. Let us pick 
an arbitrary real unit vector n and take a full convolution of riirij with the above expression: 

muj K'\k) = {\n-u\^) , (A29) 

whence, and by using (A6), we obtain, for all k, 

< ninj«;'^ (k) < (|u|^) = K"(k), (A30) 
Cn • kl^ 

0<K{k) + k{k)^-j-^<dK{k)+K{k). (A31) 

Wc may now write (n • k)^ — fc^ cos^ 9, where 6 is the (arbitrary) angle between k and n. Thus we arrive at 
the following realizability theorem: for any value of 9 and any k the following two inequalities must hold: 

K(fc) + kik) cos^ 9>0, (A32) 
{d-l)K{k) + k{k)sm^9>0. (A33) 

In practice, we will be mostly content to use just the particular cases of the above arising from setting ^ = 
and 6 = 7r/2: 

K{k) + k{k) > 0, (A34) 
K{k) > 0. (A35) 

Let us now derive the realizability constraints that operate in the configuration space. Consider the 
two-point structure function of the field w': 

( [u' (y) - (0)] [u^ (y) - (0)] ) = 2k'^ (0) - 2k'^ (y) , ( A36) 

take the full double dot product of the above expression with an arbitrary unit vector n, and use the 
expression ( A4) for the correlation tensor k^^ . The resulting realizability theorem is 

- KLL{y) cos2 9 - siu^ ^ = 1 ( |n • [«'(y) - u\0)] |') > (A37) 

for any y and any 9 (the angle between n and y) . The inequalities of most practical value are again obtained 
for 6* = and 9 = 7r/2: 

KLLiv) < kll(O), (A38) 

KNN{y) <{d- l)AtLL(O) = KJVW(O). (A39) 

Thus, the two-point correlation functions can never exceed the values they take when the two points coincide. 
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A.4. Small-Scale Expansion of Second-Order Correlation Functions 

Since much of this work is concerned with one-point statistics, let us establish a correspondence between 
one- and two-point correlation properties of our field u*. For y — > 0, we will routinely expand 

i^'Hy) = i^oS'' - i «2 {y^S'^ + 2afy') + i n^y^ {yH'^ + 2hy'y^) +■■■ . (A40) 

Here a and b are the so-called compressibility parameters which change within the intervals 

1 2 

< a < 1 and - - — - < 6 < 2, (A41) 



d+l~ ~ + 3 

where the lower bounds correspond to the incompressible case, the upper bounds to the irrotational one. 
Making use of the formula (All), we expand the expressions (A14) and (A15) for K,LL{y) and K,NN{y) at 
small y, substitute them into (A4), and, by comparing with (A40), establish the small-scale expansions of 
the velocity correlation functions in the x space, 

HLL{y) = Ko- ^(l + 2a)K22/' + i(l + 26)Av4/ + --- , (A42) 

KNwiy) = {d-1) (^Ko - ^ K2y'^ + ^ K4y'^ H ^ , (A43) 

and the expressions for the parameters of the correlator (A40) in terms of the spectral characteristics of the 
velocity field, 

1 f Adu 

" = '"- d(. + 2Kd + 4) /(S'^'"'^'- 

where the k-space integrations of radial functions could, of course, have been written more explicitly as 



{2nY (27r)'' Jq 



dkk'^-\ (A49) 



Finally, let us remark that the expression (14) that was obtained in Sec. 2 for the initial growth rate of 
the total magnetic energy can be, in view of the above, written as follows: 

7 = ^ (rf - 1) (da + 2)k2 = ^ 1^ (d + 2 + /3) K2, (A50) 

where /3 = d[l + a{d+ 1)]. The above expression is naturally in agreement with the results known in the 
one-point setting (see, e.g., Boldyrev and Schekochihin 2000; Schekochihin and Kulsrud 2001). Another 
convenient expression (or definition) for 7, that we use in Sec. 3, is 

7=^(l + 2a).2 = ^%P, (A51) 
where A is defined by the formula (22). 



-35- 



A.5. Angle-Integration Identities 

Dealing with d-dimensional wave- number integrations while deriving the mode-coupling equation (11), 
the evolution laws (13) and (16), and the SSF equation (21), required repeated use of some simple angle- 
integration identities which we list here for completeness: 



// 



/ 



/ 



/ 



dCldniTij = (A53) 

dfldriinjnkni = 2) ('^^J'^'^' ^ik^ji + ^n^kj) , (A54) 

dQddn'ain-n'f = (A55) 

j jdUddU'M-r^f = (A56) 

where dfid = A^%va.Q\Adi . . . sm'^~'^ 9d~2d0d-2 is the d-dimensional angular differential, and n = k/A; is a 

unit vector in the direction of k. The above identities can be derived in various ways, one of them being via 
formula (A9). 



B. THE FURUTSU-NOVIKOV FORMULA 

Here we state an extremely useful Gaussian-averaging theorem that is widely used throughout this work. 
This result is due to Funitsu (1963) and Novikov (1965) and is generally referred to as the Furutsu-Novikov, 

or "Gaussian-integration," formula. 

Consider a random Gaussian vector field u\q), where q is the aggregate of all variables it depends on 
(such as time, spatial position, Fourier variable, or vector index). Let R[u] be some functional of the field u*. 
Then the following holds: 

{u\q)R[n]) = J dq' {u\q)u\q')) (^f|^) , (Bl) 

where the integration is over all possible values of q'. 

This property of Gaussian fields is simply a generalization of the well-known "splitting" property of the 
Gaussian averages. 



C. LARGE-SCALE CORRECTIONS TO THE INERTIAL-RANGE ASYMPTOTICS 

In this Appendix, we work out the first "large-scale" corrections to the inertial-range asymptotics (77) 
and (78) of the Kazantsev particle's mass and the potential it lives in, as well as the appropriate corrected 
asymptotic forms of the solutions. These corrections arc needed in order for us to get a better idea of 
the general form of the Kazantsev quantum mechanics. They were also appealed to when the asymptotic 
matching was performed in § 3.3. 
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In order to obtain the desired corrections, wc mnst use the expansion (A40) of the velocity correlator (50) 
up to the fourth order in y. The fourth-order Taylor expansions of the renormalized diffusivities are 



K{y) 
Q{y) 



2r? + ^ (1 + 2a)K2t/^ - ^ (1 + 26)'^4y^ 



2?? + ^ K22/^ 



1 4 



and we find the following corrected asymptotics for mass and potential in the inertial range: 

1 



V{y) ~ ~Vo + V2v\ 



1 1 + 26 K4 , 
1 + r z 



2 1 + 2a K2 



(CI) 

(C2) 

(C3) 

(C4) 



where V2 = - 



6(l + 2a)2 



K4,. 



The correction V2y^ in the formula (C4) is positive in the 3D incompressible case, changes sign at 6 = 1/2, 
and remains negative for all b> 1/2, including in the irrotational case. This means that the value V{y) ~ —Vq 
represents an intermediate "plateau," and not the bottom of the potential well, which extends further down 
[how far down is, of course, impossible to determine on the basis of the Taylor expansion (C4); a representative 
example of such a potential is given in Fig. 9] . The danger thus arises that the Kazantsev particle might 
slide into this deeper recess and thereby destroy the small-scale universality, upon which hinged the validity 
of the results obtained in § 3.3. Whether or not this happens in reality may depend on the large-scale 
structure of the potential and on the interplay between it and the structure of the variable mass m(y). The 
importance of the latter effect is clear from the form of the eigenvalue problem (64) itself, which can be 
reinterpreted as a Schrodinger equation for a particle with constant mass and zero energy in a A-dependent 
effective potential Ux{y): 



-< + i7A(y)V = 0, 
Ux{y) = 2m{y)[V{y) + \^]. 



(C5) 
(C6) 



We must now treat A as a parameter to be chosen in such way that the zero-energy eigenstate exist and 
be a ground state. The effective potential U\{y) explicitly incorporates both influences [variable mass and 
potential V{y)] that act on the Kazantsev particle. In the inertial range, it is, of course, an inverse-square 
potential. Due to the inverse-square decay of m(y), the fourth-order effects lead only to a constant correction, 
namely. 



Ux{y) 



A 



■C/2, 



(C7) 



where U2 = 



1 



V2--{1 + 2b)K4 J 



The solutions of Eq. (C5) with the effective potential (C7) are the modified Bessel functions: 
These are the solutions we referred to in justifying the asymptotic (91) of § 3.3. 



(C8) 



The potentials V{y) and Ux{y) for a particular 3D irrotational velocity field with the transverse correla- 
tion function chosen to be KjvjvCy) = ^Kq exp (— «;22/^/2) , are shown in Fig. 9 and Fig. 10. In this specific case. 
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the straightforward numerical solution of the eigenvalue problem (64) proved that the Kazantsev particle was 
"heavy enough" to settle on the flat "porch" corresponding to the inertial-range asymptotic V{y) ~ — Vq, 
rather than slide off towards larger scales and deeper into the potential well. The small-scale universality 
therefore persisted. 
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Fig. 1. — Sketch of the smaU-scale dynamo in 2D. The number of excited modes grows, while each individual 
mode decays. Both processes occur exponentially fast. The total growth rate of the magnetic energy in 
the diffusion-free regime is 27. The per-mode decay rate and the spectral index depend on the degree of 
compressibility. The numbers on the sketch correspond to the incompressible case. 
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Fig. 2. — Sketch of the small-scale dynamo in 3D. Both the number of modes and the amplitudes of all 
individual modes grow exponentially in time. The total growth rate of the magnetic energy in the diffusion- 
free regime is 2j. The per-modc growth rate depends on the degree of compressibility, but the spectrum is 
always k^^"^. The numbers on the sketch correspond to the incompressible case. 
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Fig. 3. — Dependence of the growth rate of each individual mode, Ao, on the compressibility parameter a 
in 2 and 3 dimensions [see formula (28)] and Table 2. 
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Fig. 4. — The log-log plot of the spectrum (35) in 3D incompressible flow. The spectrum M{t, k) is nor- 
malized to the total energy W{t) and shown at times 37~^, 5j^^, 77^^, and 97^^. For reference, the line 
corresponding to the expected k^^^ inertial-range spectrum is also plotted. The initial wave number was 
taken to be ko = 0.05k„. All wave numbers are normalized to the resistive wave number k„- 
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Fig. 6. — Evolution of the mean square wave number of the magnetic fluctuations, [fc^(t)] , in two and 
three dimensions in the case of incompressible flow. The time is measured in imits of [defined by the 
formula (14)]. The wave numbers are measured in units of the inverse resistive scale /c^. The initial wave 
number was taken to be ko = O.OSfcn. 
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Fig. 7. — A schematic plot of the variable mass in Kazantsev quantum mechanics. 
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Fig. 8. — A schematic plot of the potential in Kazantsev quantum mechanics in three dimensions [see 
formula (62)]. In its particulars, this plot depicts the potential V{y) for the 3D incompressible flow. 
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Fig. 9. — The Kazantsev potential in the case of 3D irrotational flow. For reference, the Hne corresponding 
to V = — Vo is also shown. V{y) is measured in units of Vb, y in units of {^0/^2)^^^ ■ The transverse 
velocity correlation function was chosen in the form KNN{y) = 2Koexp (— K22/^/2). Due to the irrotational 
property of the flow, the longitudinal correlation function had to be calculated according to the formula (A21) 
of Appendix A. 2. 
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Fig. 10.- The form of the function U\{y) — 2m{y)\y{y) + A7] which corresponds to the potential V{y) 
plotted in Fig. 9 and A = Aq = 1/8. The bottom of the potential well extends to U\ ~ —k^/t] and is not 
shown. 
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Table 1. Physical parameters in the Galaxy and Protogalaxy. 



Parameter 


Notation 


Galaxy (warm ISM) 


Protogalaxy 




Densities and Temperatures 


neutral (hydrogen) density, cm~ 
ion (proton) density, cm~^ 
temperature, K 


^ «„ 

rti 

T T-i T'g 


1 
1 

10* 




10"^ 
10<^ 


Length Scales, cm (3 - 10^^ cm ?s 1 parsec) 


system size 

magnetic field coherence 

energy-containing 

viscous cutoff 

neutral mean free path 

ion mean free path 

resistive cutoff 

ion skin depth 


L 

~ "i^th rnin{rTiTi, r„i} 
~ ^th min{ri^, Tin} 


10== 

10=1 

10=° 

lO^** 

10" 

10^= 

lO'' 

10^ 


10=* 

unknown 

10=* 

10=° 

10i« 

lO'' 

10« 


Time Scales, sec (tt ■ 10^ sec fa 1 year) 


system lifetime 
system rotation 
largest-eddy turnover 
viscous/smallest-eddy turnover 
neutral-neutral collision 
ion-neutral / neutral-ion collision 
ion-ion collision 


T -1 

Teddy (fcixModdy) 

T„„ ~ 47rr=jjjj^t)thn7i 


IQi^ 

IQi^ 

10" 

10^= 

10^ 

10"* 

10'^ 


10^^ 

10^" 
10" 

10^= 


Velocities, cm/sec 


thermal / shock/largest-eddy 
smallest-eddy 


-"th - {kBT/rrn)^^^ 
"eddy ~ Re~^^'^Vth 


10^ 
10^ 


lO'' 
10** 


Viscosities, cm^/sec 


neutral viscosity 
ion viscosity 

magnetic (Spitzer) diffusivity 


Vr^ ~ t;?hT„„ 
I^i ~ fthT,» 

77 ~ {kBT)-^''^ml'^e^c'^ In A/4;r 


10=1 
10^^ 
10^ 


10=s 
10* 


Dimensionless Numbers 


hydrodynamic Reynolds 
magnetic Reynolds 
magnetic Prandtl 


Re vtbL/v 

Rm ~ va^L/-q 

Pr ~ u/ri ~ (kn/k^ f 


10^ 

10i« 

10" 


10* 

10='' 

10== 
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Table 2. Parameters of the SSF equation. 



Dimension Velocity Field a A B C sq Xq X, 



d = 3 Incompressible -1/4 1/5 2/5 4/5 -1/2 3/2 3/4 4/5 

Irrotational 1 3/10 3/5 1/5 -1/2 3/2 1/8 1/5 

d = 2 Incompressible -1/3 1/4 -1/4 1 2 -1/4 -3/16 

Irrotational 1 3/4 5/4 -1/3 2/3 -1/12 5/48 



